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COMPUTATIONAL SOLUTION OF ATMOSPHERIC CHEMISTRY PROBLEMS 


By 

Jawed A. Jafri 1 and R. L. Ake 2 
INTRODUCTION 

Extensive studies were carried out on problems of interest in 
atmospheric chenistry during the grant period. In addition to several minor 
projects, four major projects were carried out during the course of this 
research and are described in the following sections. 

I. THEORETICAL STUDIES OF GROUND AND LOW-LYING EXCITED STATES OF Cl U 2 

Chlorine introduced into the atmosphere as a result of natural or an- 
thropogenic processes reacts with ozone via 

Cl + o 3 ► CIO + 0 2 

Cl + 0 2 - Cl 00 

and 

CIO + CIO Cl 00 + Cl 

CIO + 0 2 + M— - C10 3 + M 

Detailed computations were carried out on the C10 2 molecule. Ground 
and several low-lying excited-state potential energy curves for the doublet 
and quartet states in A' and A" symmetries were calculated. The potential 
energy curves and other properties computed for various states were used 

*NASA Langley Research Center, Hampton, Virginia 23665 
Associate Professor, Department of Chemical Sciences, Old Dominion 
University, Norfolk, Virginia 23508. 


to determine the nature of bonding between the Cl and 0 2 fragments, the 
wavelengths of the solar spectrum which would cause photolysis, and the 
photolysis products. The results of this investigation were published in 
the Journal of Chemical Physics . A reprint of the article detailing and 
describing the results is included in the Appendix. 

Investigations were carried out on the other product, Cl 0 3 . Double 
zeta or better quality basis sets with polarization functions were employed. 
C10 3 with this basis became too large to handle with the software programs 
available at that time. SCF calculations were carried out at several 
geometries to find the equilibrium geometry, but no binding was obtained at 
this level. Some GVB calculations were also carried out, but these 
calculations were very large and also expansive. In order to study this 
system properly, large scale MCSCF and Cl calculations are needed together 
with several parameter optimizations, and since these needs exceeded the 
resources at that time, these calculations were postponed. 
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II. GROUND AND EXCITED STATE POTENTIAL ENERGY SURFACES OF THE 
METHYL PEROXY RADICAL, CH 3 0 2 

The methyl peroxy, CH 3 0 2 , radical is an important intermediate in 
combustion processes and in the degradation of methane in the atmosphere. 
CH 3 0 2 and subsequent products of methane degradation participate in chain 
reactions which are important in the production and destruction of ozone in 
the troposphere and stratosphere. The fractional abundance of CH 3 0 2 in 
the atmosphere is about 10 - 11 with a lifetime of 10 3 seconds (Refs. 1, 3). 
Its role in the chain reaction involving methane with HO, NO, and 0 2 
molecules might indicate an increase in ozone concentration near the tropo- 
sphere where the NO concentration may be significant (Ref. 2). The prod- 
uct of self reaction of CH 3 O 2 radicals subsequently reacting with 0 2 
yields HCHO, H0 2 and CH 3 0 radicals and is non terminating. The aosorp- 

tion spectrum of CH 30 2 in the gas phase has been studied extensively. The 
only theoretical ab initio work is due to Bair and Goddard (Ref. 4), who 

examined the ground state of CH 3 O 2 to determine the equilibrium geometry 
using a valence double zeta basis set augmented by polarization functions in 
a GVB-CI scheme (Ref. 5). These computations, however, do not determine the 
dissociation energy or the nature of the excited states. 

Work, under this grant, was carried out to calculate the ground state 
geometry, the dissociation energies, the nature of excited states and photo- 
lysis products using a larger basis set and a more extensive Cl scheme. Two 
sets of computations were carried out which employed the same MCSCF proce- 
dure but different basis sets and different Cl schemes. In the first set of 
calculations a double zeta quality Cartesian Guassian basis set was chosen 
consisting of Dunning's (Ref. 6) [4s, 2p/2s] contraction of Huzinaga's (Ref. 
1) (9s,5p/4s) components, with a scale factor x, - 1.2 chosen for 
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hydrogen's functions. The SCF solutions were obtained for the single and 
multiconfiguration wave functions using a modified version of the GVB pro- 
grams due to Goddard and co-workers (Ref. 8). Double excitations to all 
antibonding orbitals from the bonding orbitals, including C-H bonds, in 
CH 3 0 2 , were allowed in the MCSCF procedure, and the set of orbitals thus 
obtained was used in this first set of calculations. These "pseudo first 
order" Cl calculations utilized CSF selection, by the A k (Ref. 9) method, 
from a list generated by considering all valence double excitations from a 
multiple CSF reference function with the limitation that only one electron 
could occupy a virtual orbital. These calculations were used to carry out a 
preliminary search for the lowest three states in 2 A" and 2 A' symne- 
tries. 

The second set of computations on CH 3 0 2 utilized the sane MCSCF pro- 
cedure as above with an augmented basis set. Polarization was provided by a 
set of d-type functions on carbon (a = 0.75), and oxygen (a = 0.85) and 
by a set of p-type functions (a = 1.0) on the hydrogen atoms. A pseu-do 
first order technique, similar to the one in the first set of calculations, 
was used to obtain a set of approximate natural orbitals. These ANO's 
were subsequently used to carry out a selected singles and doubles Cl on 
the ground state. In this set of calculations, CSF's were selected from a 

set generated by single and double excitations into the valence and virtual 
space from a multiple CSF reference function, using the AK method. 

The potential energies for the three lowest 2 A" and 2 A' states for 

various CH3O2 geometries, where the R(0-0) coordinate was varied keeping 
the R(C-O) at 2.725 bohr, are presented in Table 1.* These energies are 
plotted in Fig. 1.* The energies and dipole moments of the lowest 2 A" and 

*Tables and figures cited in each section appear at the end of the section. 
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2 A‘ state as a function of R(C0) [R(0-0) = 2.783 bohr] are presented in 

Table 2 and plotted in Fig. 2. These calculations utilized a double zeta 
basis set and "pseudo first order" Cl. The first excited state in 2 A" 
symmetry along the C-0 bond stretch coordinate (Fig. 2) which should domin- 
ate the photochemistry of CH 3 0 2 , exhibits a bound region near R(C-O) = 
2.75 bohr and an avoided crossing with the third 2 A" state near R(C-O) = 

3.5 bohr. Along the 0-0 stretching coordinate this 2 A" state is essen- 
tially repulsive near R(0-0) = 2.53, but exhibits slight bonding at longer 
bond lengths. The shape of the 2 A" surface indicates that CH 3 O0 should 
be the principle photolysis products. If this state is slightly bound at 
long 0-0 bond lengths, a temperature dependent quantum yield for photoly- 
sis could be expected in the long wavelength tail of the absorption band for 
this state. 

The first 2 A' state is bound in both the 0-0 and C-0 bond stretching 
coordinates as expected on the basis of spectroscopic investigation of 
CH 3 0 2 absorption in the near infrared. The second 2A ' state is r e p U lsive 
in both bond stretching coordinates and exhibits a crossing with the second 
2 A" state on the outer limbs of this state's well on the 0-0 bond 
stretching coordinate. The location of the crossing on this particular 
slice through the surface indicates that the resulting predissociation of 

2 2 A", CH 3 O 2 to CH 3 O+O should not be important. This indication, 
however, is not definitive due to the sparse coverage of the region and the 

approximations in the level of Cl treatment. The excitation energy of the 
2 2 A' state in the region of the ground state equilibrium should prevent it 

from being significant in the atmospheric photochemistry of CH 3 O 2 . 

The second series of computations with more comprehensive Cl treat- 
ment and larger basis set as discussed above was carried out. The energies 
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for the first two 2 A" states along the R(0-0) coordinate with R(C-O) 
fixed at 2.7483 bohr are presented in Table 3 and plotted in Fig. 3. The 
energy of the second 2 A" is lower (at long bond lengths) along the R(0-0) 
stretching coordinate and exhibits bonding similar to that found in the 
"pseudo first order" calculation along that slice in this coordinate. The 

dissociation energy for CH 3 O 2 * CHO 3 + 0 determined at this level of 
calculation is 3.36 eV. 

The refined computations on the potential energy curves along the 
R(C-O) stretching coordinate, with R(0-0) fixed at 2.5606 bohr are pre- 
sented in Table 4 and plotted in Fig. 4. The computations give similar 
results as in "pseudo first order' calculations with the avoided curve 
crossing between second and third 2 A" states near R(C-O) = 3.2 bohr. The 
calculation reported in Tables 3 and 4 utilized different approximte natural 
orbitals for the first and second roots. The dissociation energy for 
CH 3 0 2 -»- CH 3 + 0 2 at this level of computation was 2.08 eV. 

The potential energy curve along a 45° angle in the plane defined by 
the R(C-O) and R(0-0) coordinates is presented in Table 5 and plotted in 
Fig. 5. Here again, different approximate natural orbitals were employed 
for the first and second root calculations. The second 2 A" state is 
purely repulsive along this coordinate. 

The vertical excitation energy between the x 2 A" «- 1 2 A' is computed 
to be 1.23 eV in the "pseudo first order" type calculations. 
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Table 1. Energies of the lowest 2 A" and 2 A' states of CH 3 0 2 as a 
function of R(0-0); R(C-O) fixed at 2.725 bohr. 


R(0-0) 

State 

El 

e 2 

e 3 

2.1505 

2 A" 

-189.298021 

-188.984950 

-188.927101 


2 A* 

-189.239457 

-188.974894 

-188.968584 

2.430 

2 A" 

-189.366776 

-189.137986 

-189.006940 


2 A' 

-189.323150 

-189.108858 

-189.051324 

2.530 

2 A" 

-189.375718 

-189.170388 

-189.060038 


2 A' 

-189.336555 

-189.137441 

-189.102755 

2.680 

2 A" 

-189.380272 

-189.209729 

-189.121958 


2 A' 

-189.347146 

-189.167850 

-189.161750 

2.783 

2 A" 

-189.379000 

-189.227709 

-189.153899 


2 A' 

-189.349894 

-189.191495 

-189.182926 

2.900 

2 A" 

-189.375016 

-189.242442 

-189.181881 


2 A' 

-189.349249 

-189.217065 

-189.194490 

3.036 

2 A" 

-189.368279 

-189.253157 

-189.205724 


2 A' 

-189.347488 

-189.239335 

-189.203314 

3.289 

2 A" 

-189.353644 

-189.260269 

-189.230523 


2 A' 

-189.339151 

-189.264805 

-189.212428 

3.619 

2 A" 

-189.336507 

-189.267475 

-189.251220 


2 A 1 

-189.327569 

-189.281986 

-189.217849 

3.796 

2 A" 

-189.329485 

-189.264481 

-189.252198 


2 A' 

-189.324086 

-189.288061 

-189.216604 

00 

2 A" 

-189.309762 

-189.295345 

-189.224714 


2 A' 

-189.304290 

-189.302305 

-189.216360 
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as a 


Table 2. 

Energies of 
function of 

the lowest 2 A" and 2 A' states of 
R(C-O) ; R(0-0) fixed at 2.783 bohr. 

CH 3 0 2 as a 

R(C-O) 

State 

Ei 

e 2 

e 3 

2.4425 

2 A" 

-189.350966 

-189.213618 

-189.120739 


2 A‘ 

-189.321927 

-189.165722 

-189.161895 

2.60 

2 A" 

-189.371913 

-189.226331 

-189.151162 


2 A' 

-189.342519 

-189.183436 

-189.179728 

2.725 

2 A" 

-189.379219 

-189.228451 

-189.158422 


2 A' 

-189.349362 

-189.190420 

-189.183653 

2.85 

2 A" 

-189.380861 

-189.226953 

-189.160862 


2 A‘ 

-189.350579 

-189.191848 

-189.183186 

2.998 

2 A" 

-189.378283 

-189.221670 

-189.159368 


2 A‘ 

-189.347174 

-189.188517 

-189.178544 

3.270 

2 A" 

-189.365512 

-189.207281 

-189.154780 


2 A' 

-189.332228 

-189.175355 

-189.164866 

3.543 

2 A" 

-189.347680 

-189.191854 

-189.168744 


2 A' 

-189.312382 

-189.156824 

-189.147050 

3.815 

2 A" 

-189.329582 

-189.200697 

-189.172280 


2 A' 

-189.294127 

-189.164362 

-189.130020 

oo 

2 A" 

-189.315773 

-189.089019 

-189.981690 
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Table 3. Energies (in Hartrees) of the first two 2 A" states along the 
R(0-0) stretching coordinate with R(C-O) = 2.7483 bohr. 


R (0-0) 

Ei* 

e 2 ** 

2.411 

-189.627009 

-189.360728 

2.5606 

-189.634754 

-189.399470 

2.711 

-189.629232 

-189.429652 

2.80 

-189.624619 

-189.438978 

2.861 

-189.619795 

-189.460298 

3.00 

-189.609292 

-189.448384 

3.30 

-189.583719 

-189.476526 

3.80 

-189.554940 

-189.444410 

00 

-189.511096 



*Energies calculated using approximate natural orbitals for the 2 A" state. 
**Energies calculated using approximate natural orbitals for the second 2 A" 
state. 
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Table 4. Energies (in Hartrees) , for the first two 2 A" states in the R(C- 
0) stretching coordinate. R(0-0) was fixed at 2.5606 bohr. 


R(C-O) 

El* 

e 2 ** 

2.598 

-189.631742 

-189.408465 

2.748 

-189.634754 

-189.399470 

2.898 

-189.629929 

-189.395334 

3.10 

-189.615957 

-189.380685 

3.30 

-189.587609 

-189.363211 

3.60 

-189.575821 

-189.402720 

3.90 

-189.556006 

-189.438198 


* Energies calculated using approximate natural orbitals for X 2 A" 
** Energies calculated using approximate natural orbitals for 2 2 A" 


state. 

state. 
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Table 5. Energies (in Hartrees) of the X 2 A" and 2 2 A" states along the 
stretching coordinate along a 45° angle through the plane defined 
by R(C-O) and R(0-0) stretches. The first and second root 2 A" 
states employed approximate natural orbitals determined 
separately. 


R 

Ei 

e 2 

3.5443 

-189.627012 

-189.363965 

3.6033 

-189.632922 

-189.380388 

3.7583 

-189.634754 

-189.399470 

3.9133 

-189.632372 

-189.415834 

3.9683 

-189.626890 

-189.422632 

4.1803 

-189.606836 

-189.425222 
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189.121 










Energy (Hartrees) 


189.35 


189.40 


189.45 


189.50 


-189.55 


-189.60 


-189.65 



2.4 2.8 3.2 3.6 4.0 

R(0-0) in bohr 

Figure 3. Potential energy of the X 2 A" and 2 2 A" states of CH 3 O 2 
along the R(O-O) coordinate. The calculations used a 
double zeta basis augmented by polarization functions and a 
selected single and double Cl in approximate natural 
orbital basis. R(C-O) fixed at 2.7483 bohr. 
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189.30 


•189. 35 J 
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•189.50 4 


189.55 


-189.60 



-189.65 
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Figure 5. 


3.756 3.968 4.180 

R (bohr) 

Potential energy of X 2 A" and 2 2 A" states of CH 3 O 2 
along 45° angle in the plane defined by R(C-O) and 
R(0-0) coordinates. The calculations used a double 
zeta basis augmented by polarization functions and selected 
singles and doubles Cl in approximate natural orbital 
basis. 


III. ELECTRONIC STATES OF THE FO RADICAL. 

The ground and lower excited states of the FO radical have been 
investigated using a large basis set with MCSF and Cl techniques. FO has 
long been an elusive molecule. Direct evidence of its existence was not 
obtained until 1964 when Arkell, Reinhard and Larson (Ref. 10) detected the 
IR band at 1028 cm- 1 in an argon matrix containing the photolysis products 
of F 2 0. Andrews (Ref. 11) subsequently studied the fundamental bands using 

Raman spectroscopy of F 2 O photolysis products in an Argon matrix. Early 
estimates of dissociation energy of FO ranged from 24 Kcal/mole based on 
direct evidence from electron impact studies to 46 Kcal/mole inferred from 
the thermolysis of F 2 0. Later experiments produced values for the dis- 
sociation energy of FO as high as 59 Kcal/mole. 

O'Hare and Wahl (Ref. 12) and Lathan et. al (Ref. 13) have reported 
Hartree Fock calculations on the ground state of the FO radical. 

Langhoff, et. al (Ref. 16) studied the X 2 w potentials and dipole moment 
functions for FO at the CASSCF, externally contracted Cl and MRCI(SD) 
levels using double zeta plus polarization and extended Gaussian basis sets. 
No experimental or theoretical results are available for the excited states 
of FO radical. 

Ab initio calculations using an accurate basis and including electron 
correlation, on the ground and excited of FO radical were carried out. 
MCSCF and first order Cl calculations have been carried out for the ground 
and excited states, and selected singles and doubles Cl for the ground state 
obtaining potential energy curve, dipole moment curves, and transition mo- 
ments between various excited states and X 2 ir ground state. The inten- 
sities of the vibrational transitions have been obtained in order to provide 
a basis for infrared measurements of the concentration of FO in either the 
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laboratory or the upper atmosphere where the radical may be produced by the 
photochemistry of flurocarbons. 

The basis set used in the ground state and all excited state calcula- 
tions was a triple zeta contraction of Van Dujneveldt's (Refs. 14, 15) 
(llS,7p) valence augmented by two sets of d functions. The exponents for 
the d-type functions on the oxygen atom were 2.091 and 0.6194 and on the 
fluorine atom were 2.214 and 0.656. This basis set yields an SCF energy for 
the X 2 tt state at 2.5 bohr internuclear distance which is 3 x 10- 4 Hartrees 
higher than O'Hare and Wahl's near Hartree-Fock value. 

The initial orbitals were obtained by carrying out a GVB calculation 
where the excitation involved replacement of the FO bonding orbital with the 
antibonding orbital. The orbitals arising principally from combinations of 
the atomic Is orbitals were frozen for all subsequent calculations. Tne 
orbitals were further optimized by carrying out a full valence MCSCF calcu- 
lation using a general second order MCSCF technique developed Dy Lengsfield 
(Ref. 16). These calculations were carried out in C ^ symmetry. The vir- 
tual orbitals of A symmetry which arise from the d-orbitals were omitted 
from the full valence MCSCF calculations in order to prevent their being 
mixed with the orbitals of the w symmetry. These A type orbitals were 
included in the subsequent Cl calculation. The Cl calculations utilized 
the full symmetry of the molecule. The ground state of the FO radical was 
investigated at two levels of Cl treatment. The first order Cl treatment is 
generally regarded as providing accurate one electron properties. This was 
followed by an (Ref. 8) selected single and double excitation, multi- 
reference (SD-CI) treatment in order to determine the degree of convergence 
of the dissociation energy. To ensure a uniform treatment throughout the 
potential energy surface, important CSF's were selected from the MCSCF 
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results at the equilibriun geometry and for larger separation of the 0 and 
F atoms. The union of these two sets of important CSF's were used across 
the potential energy curve. The vibrational Schrodinger equation was inte- 
grated using the calculated potential energy curve. 

The matrix elements of the dipole moment operator were calculated in 
order to determine the vibrational ly averaged dipole moment and vibrational 
intensities. Separate MCSCF calculations were carried out for each excited 
state investigated. Full valence and first order computations were carried 
out for the full curve on 2 2 ir and 3 2 ir excited states and at the ground 
state equilibrium geometry for other doublet and quartet states. These 
results are presented in Tables 6-16, and are plotted in Figs. 6-10. 

The potential energy curve for the low lying excited states of FO 
radical were found to be repulsive, and an avoided crossing was found be- 
tween 2 2 ir and 3 2 tt states near 3.3 bohr. The repulsive nature of the 
excited states provides an explanation for the failure of experiments in- 
tended to detect band spectrum for the FO radical. The intersection De- 
tween 2 2 tt and 3 2 P states may lead to the production of small amounts 

of atomic oxygen in the state if the FO radical is in fact produced 

and photodissoci ated in the upper atmosphere. The calculated excitation 
energies (F.O. Cl) indicate that the absorption continuum due to the second 
2 ir state should have a maximum at approximately 6.4 eV. This absorption 
maximum due to the weaker transition of the 2 a and 2 i- state should 
occur at about 5.9 eV, and the weak absorption due to the 2 e + state should 
peak at about 6.6 eV. These spectra should combine to yield a relatively 
broad and featureless continuum. 

The dipole moment curve for the ground state crosses zero near 2.5 bohr 
separation and the magnitude of the dipole moment operator supports the 
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conjecture that the failure of the experimental attempts to measure the epr 
spectrum was due to a small dipole moment. The dipole moment for the 2 2 tt 
state goes to zero when the x 2 ir dipole goes through a maximum. The theo- 
retical dissociation energy and vibrational frequencies agree well with the 
available experimental measurements. 


Table 6. Nature of Ground State Wave function for FO at Various Levels of 
Computation. 

Coefficient at Coefficient at Coefficient at 

Occupation^ a ) MCSCF level F0CP b ^ level SD^-Selected Cl level 


(core) 

5cr2 

1tt4 

2tt 3 

0.9671 

0.9477 

0.9484 

(core) 

5cr2 

lir 3 

2tt4 

0.1141 

0.1599 

0.0903 

(core) 

5a 

6a 

lit 3 2ir** 

0.1189 

0.1103 

0.0798 

(core) 

6a 2 

lir*» 

2tt 3 

0.1762 

0.1403 

0.1409 


R(F-O) = 2.6 bohr 

(a) Only one spin coupling was important for each occupation. The 
occupation is used to represent the most important CSF for that 
occupation. 

(b) First Order Cl 

(c) Singles and Doubles selected Cl 
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Table 7. FO Ground State Energies (Hartrees) as a Function of Internuclear 


Separation 

at Various Levels 

of Computation. 


R(F-O) (bohr) 

MCSCF 

FOCI 

Selected S-D Cl 

2.10 

-174.180888 

-174.291148 


2.30 

-174.236725 

-174.349259 

-174.584322 

2.35 

-174.244387 

-174.357284 

-174.592071 

2.40 

-174.250186 

-174.363248 

-174.597805 

2.45 

-174.255001 

-174.370319 

-174.603071 

2.50 

-174.257808 

-174.370319 

-174.605677 

2.55 

-174.259504 

-174.371665 

-174.607192 

2.60 

-174.259872 

-174.371889 

-174.608246 

2.65 

-174.260276 

-174.371350 

-174.607299 

2.70 

-174.259651 

-174.369992 

-174.606288 

2.75 

-174.258521 

-174.368007 

-174.604676 

2.80 

-174.253093 

-174.365510 

-174.602483 

2.85 

-174.256992 

-174.362598 

-174.600048 

2.95 

-174.250874 

-174.355873 

-174.594218 

3.00 

-174.248560 

-174.352177 

-174.593418 

3.05 

-174.246202 

-174.348424 

-174.587705 

3.10 

-174.243845 

-174.344577 

-174.584321 

3.30 

-174.235069 

-174.329412 

-174.570593 

3.70 

-174.223517 

-174.305429 

-174.547950 

4.10 

-174.218899 

-174.292375 

-174.533233 

4.50 

-174.217370 

-174.286567 

-174.524961 

4.90 

-174.216870 

-174.284135 


6.50 

-174.216860 

-174.282520 

-174.515506 
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Table 8. Energies (Hartrees) for 
at MCSCF Level. 

the 2 ir Ground 

and Excited 

States of FO 

R(F-O) (bohr) 

x 2 i r 

2 2 tt 

3 2 ir 

4 2 TT 

2.1 

-174.180888 




2.3 

-174.236725 




2.5 

-174.257808 

-174.00420 



2.6 

-174.259872 

-174.03239 

-173.862516 


2.7 

-174.259651 

-174.05410 

-173.915170 


2.9 

-174.253093 

-174.0839 

-173.99278 


3.0 

-174.248560 

-174.094221 

-174.030919 


3.1 

-174.243845 

-174.102334 

-174.062809 


3.2 


-174.108807 

-174.089343 

-174.062341 

3.3 

-174.235069 

-174.11390 

-174.109419 

-174.076966 

3.4 


-174.129959 

-174.117929 

-174.08863 

3.5 


-174.14516 

-174.121464 

-174.097914 

3.6 


-174.157729 


-174.105293 

3.7 

-1714.223517 

-174.16810 

-174.128365 

-174.11162 

4.1 

-174.218899 

-174.20650 

-174.131681 


4.9 

-174.216870 

-174.21180 

-174.134974 


6.5 

-174.216860 

-174.214917 

-174.135980 
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Table 9. 


Coefficients of Important CSF's for 2 2 tt and 3 2 tt States of FO 
as a Function of Internuclear Separation. 

R(F-O) (bohr) 


Occupation 3 


3.1 

3.3 

3.4 

3.5 

SCF Config. 

2 2 7T 

-0.114 

-0.041 

-0.023 

-0.025 


3 2 tt 

-0.023 

-0.005 

0.076 

0.103 

k/J a* 

2 2 tt 

-0.023 

-0.036 

0.328 

-0.310 

A J 

3 2 tt 

-0.287 

-0.166 

0.026 

-0.040 

(*„.** J * 

2 2 tt 

0.003 

0.024 

0.796 

0.797 

J J 

3 2 tt 

0.857 

0.677 

-0.185 

-0.040 

(^x’^y) * a * 

2 2 tt 

-0.005 

0.001 

0.456 

0.465 


3 2 tt 

0.412 

0.364 

-0.116 

-0.036 

o + a* 

2 2 tt 

0.775 

0.730 

0.017 

0.010 


3 2 tt 

-0.035 

-0.437 

0.703 

0.695 

Q 

+ 

Q 

* 

2 2 7T 

0.087 

0.236 

0.009 

0.009 

ir ir* 

3 2 tt 

-0.023 

-0.173 

0.254 

0.276 

(a) 2 - (a*) 2 

2 2 tt 

0.519 

0.549 

0.010 

0.004 


3 2 tt 

-0.027 

-0.344 

0.504 

0.524 

(„)2 . („.)2 

2 2 tt 

-0.115 

-0.155 

-0.001 

0.010 

7T + TT* 

3 2 tt 

-0.080 

0.168 

-0.169 

-0.189 


a. In each case, only one spin coupling was important for each occupation. 
The occupations are therefore used to represent the dominant CSF. 



Table 10. 

First Order Cl 
FO. 

Energies (Hartrees) for 

it Excited States 

R(F-O) (bohr) x2tt 

21* 

32ir 

2.6 

-174.371889 

-174.13771766 

-174.9471872 

3.0 

-174.352177 

-174.18901197 

-174.10948514 

3.2 


-174.1970095 

-174.1625669 

3.3 

-174.305429 

-174.1991699 

-174.18238933 

3.4 


-174.2007387 

-174.1982876 

3.5 


-174.2141748 

-174.20198083 

3.6 


-174.2259253 

-174.20198083 

I 6.5 

-174.282520 

-174.2813276 

-174.1983615 
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Table 11. Full Valence Excited State Wave Functions of FO. 


State 


Occupation 


Component 


2 Z- 


2 E + 


2 a 


*»r- 


(core) 

(core) 

5a 2 6a Itt^tt 2 
5a 2 6a Itt 3 3 

major 

minor 

(core) 

(core) 

(core) 

5a Itt 4 2 * 4 
5a 2 6a 1W 
5o 2 6a In 3 2n 3 

minor 

major 

major 

(same 

occupation as 2 z+) 


5a 2 6a 

In 3 2n 3 

major 

5a 2 6a 

In 3 2n 3 

minor 

5a 2 6a 

In ^ 2n 2 

minor 


Table 12. 

MCSCF Energies 
States of FO. 

(Hartrees) for 


2 z + , and 

2 a Excited 



State 




R(F0) (bohr) <*z- 

2 A 

2Z+ 

2 E- 

x2tt 

2.5 

-174.0822 

-174.0043 

-173.9901 

-174.9963 

-174.2578 

2.6 

-174.1172 

-174.0423 

-174.0244 

-174.0361 

-174.2599 

2.6* 

-174.2406 

-174.1597 

-174.1302 

-174.1532 

-174.3719 

2.7 

-174.1434 

-174.0719 

-174.0556 

-174.0700 

-174.2596 

2.9 

-174.1770 

-174.1139 

-174.1015 

-174.1166 

-174.2530 

3.3 

-174.2043 

-174.1631 

-174.1583 

-174.1688 

-174.2350 

3.7 

-174.2124 

-174.1908 


-174.1953 

-174.2235 

4.1 

-174.2154 

-174.2048 



-174.2188 

4.5 

-174.2171 

-174.2112 



-174.2173 

4.9 

-174.2181 

-174.2140 



-174.2168 


♦Energies from first order Cl calculations. 
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Table 13. 


Expectation Values of the Dipole Moment 
Type b v 

A 0 

1 
2 
3 

B 0 

1 
2 
3 


a of 
<M> 
.048 
.065 
.072 
.094 
.015 
.017 
.025 
.033 


FO. 


a. In atomic units (lau = 2.5416 Debye) 

b. Calculation A used the potential and dipole moment curve from the 
selected SD-CI treatment. Calculation B used first order Cl values for 
both the potential and the dipole moment. 
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Table 14. Dipole Moments of X 2 ir 
Valence Cl Level . 

R(F-O) (bohr) 

2.5 

2.6 
2.9 
3.05 
3.30 
3.4 
3.50 
3.70 
4.1 


and the 2 tt States of F0 at Full 


X 2* 

2 2 tt 

-.0245 

-.9056 

.0418 

-.7170 

.1234 

-.7853 

.1430 

-.6387 

.1386 

-.3477 


-0.531 


-0.0484 

0.0880 

-.0397 

0.0461 

-.0249 
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Table 15. SD-CI Viorational Frequencies and Transitions i v loments a . 


v 

aG^j-i (cm-i) 

T( Av = 1) 

0 

0 


1 

941(1044) 

.012 (.013) 

2 

930(1035) 

.022 (.025) 

3 

904(1014) 

.033 (.036) 

4 

885(991) 

.040 (.047) 


a. Calculated using the potential obtained from the SD-CI. The intensities 
are in atomic units (lau = 2.5416 x 10" 18 esu cm* 1 ). The values given 
in parenthesis were obtained using the dipole, moment and potential 
energy curves from the first order calculations. The Av = 1 transition 
moments are listed with the upper vibrational state for each transition. 
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Table 16. Transition Moments of FO at Full Valance Level of Computation 
at 2.60 Bohr Internuclear Separation. 

2Z- -*• x2tt T = -0.082115 

2A ♦ x2tt T = + 0.027329 

2E+ + x2ir T = 0.17517 
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Figure 6. Potential energy curves for the x 2 ir state of FO 


06 *£ 



Figure 7. Full valence MCSCF Potential Energy Curves for Ground and Excited 
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Figure 9. Full Valence MCSCF Potential Energy Curves for Ground and Excited States 



Figure 20. x 2 ir Dipole Moment. 


IV. THEORETICAL STUDIES OF S0 2 *(H 2 0) n COMPLEXES 

The oxidation of S0 2 is an important chemical process. The phenom- 
enon of "acid rain" involves the oxidation of sulfur containing species 
including S0 2 to form SO 3 . Hydration is involved to give H 2 S0 4 which in 
turn is ultimately precipitated as sulfates (Ref. 17). The influence which 
water plays on the kinetics of oxidation is less clear. However, Rees noted 
a 4.7 x 10 3 fold increase in the formation of H 2 S0 4 from S0 2 in the 283% 
relative humidity of a smog chamber over that in the absence of water 
(Ref. 18). S0 2 may also act as the nucleation center for aerosol forma- 
tion. The exact mechanism of such a process should it exist is not known. 

The stability of a complex of S0 2 and H 2 0 was first suggested by 
Phillips from results of semiempirical calculations (Ref. 19). Holland and 
Castleman using CNDO/2 calculations suggested that an S0 2 *H 2 0 adduct 
should be stable in the gas phase and reported a binding energy of 145 mH 
and a large dipole moment (Ref. 20). 

Evidence for the existence of H 2 0*S0 2 produced from the photolysis of 
H 2 S in solid 0 2 at 15K by ultraviolet light is presented by Tso and Lee 
who conducted an FTIR study of the photolysis mixture (Ref. 21). Identifi- 

cation was based on shifted vibrational bands of S0 2 and H 2 0. 

The stability of S0 2 -water complexes has also been suggested from 
molecular beam measurements, and electric deflection experiments indicate a 
large dipole moment for the complex (Ref. 22). 

Although there have been no ab initio calculations of the S0 2 . H 2 0 

complex there are calculations for related complexes. Ab initio calcula- 
tions on S0 2 *NH 3 (Ref. 23) and S0 2 *HF (Ref. 24) show stable complexes for 
a variety of geometries. 
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An important pair of calculations were done on the closely related 

complex H 2 0- SO 3 and its rearrangement to H 2 SO 4 . A semiempirical CNDO/2 
calculation by Holland and Castleman (Ref. 25) has been followed by a recent 
ab initio calculation by Chen and Plummer (Ref. 26). 

In an effort to shed some light on a possible role of H 2 0 in the oxi- 
dation of SO 2 and elucidate the structure and energetics of formation of 
water complexes of S0 2 , we have undertaken an ab initio study. 

The study began with a reasonably sophisticated SCF-CI treatment of two 
geometries of the H 2 0.$0 2 complex. These are I and II in Fig. 11. The 
geometries of the fragments were frozen at the minimum energy theoretical 
geometries. The values used are: H 2 0 (R QH = 0.956A*, < HQH = 105.2“ ) and 

SO 2 (R so = 1.423A*, <qsq = 118.4*). The generally contracted (Ref. 27) 

Dasis sets used were: H 2 0 (9s5p/5s) /<3s2p/2s> and S0 2 (12s9p/9s5p ) !i 

<4s3p/3s2p>. Van Duijneveldt exponents (Ref. 14) were used for H and 0 
oroitals while the exponents for S were from Huzinaga (Ref. 28). 

The potential energy surface was studied at the SCF level. A configur- 
ation interaction (Cl) treatment was performed at the energy minimum of each 
conformation. Two different frozen cores were used and all single and 
double excitations outside the frozen core virtual orbitals were included. 

In CI1 the frozen core included the 2s orbital of oxygen and the 3s 
orbital of sulfur as well as the inner shell orbitals of sulfur and oxygen. 

This restriction generated 2653 configurations and 4453 cfs. In CI2 the 
frozen core included just the inner shell orbitals of sulfur and oxygen and 
generated 5547 configurations and 9735 cfs. The selection criterion in both 
cases was a cumulative discard of terms less than 5 x 10- 5 a.u. 
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The SCF potential energy surface showed a gradual slow decrease in the 
total energy as the interfragment distance was reduced, reaching a shallow 
minimum for tne two conformations, then rising steeply to become repulsive 
at distances only 0.6A° less than that which gave the minimum total energy. 
The most stable distance of conformation I was 6.1 a.u. while for conforma- 
tion II the energy minimum was reached at R = 5.5 au and r = 1.6 au with the 
S(S0 2 ) to 0(H 2 0) distance equal to 5.73 au. 

The binding energies of conformations I and II are shown in Table 17 
for the two Cl approaches. Conformation II is clearly bound for all ap- 
proaches while conformation I becomes unbound for the CI2 approach which 
substantially reduces the binding energy for both conformations over that of 
CI1. 

The results of a Mul liken population analysis (Ref. 29) are shown in 
Table 18. The SCF results for conformation II show a large movement of 
charge relative to the SCF results of the separated fragments while the 
corresponding results for conformation I show no large changes. The CI2 
approach reduced the intermolecul ar charge transfer of conformation II to 
0.001 electron. The intramolecular charge movement is also very small. 

Two basis set concerns were addressed as a part of these early calcula- 
tions. Polarization functions (3d on sulfur, 2p on hydrogen) were added to 
the basis set. The augmented basis set decreased the binding energy by 2mH 
for the particular test geometry chosen. The effect of basis set superposi- 
tion at the SCF level was found by using ghost orbitals in place of one of 
the fragments positioned at the minimum energy interfragment distance. The 
total energy was compared with a similar calculation for R = 100 au. The 
energy difference was 1.2 mH indicating only a small superposition error at 
the SCF level with the basis set used. 
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To explore the possibility of other metastable conformations of 
H 2 0*S0 2 an SCF energy gradient program (Ref. 30) was used to locate 
stable structures for the complexes. The 4-31G basis set available in the 
gradient package was used for the potential energy searches. 

Five distinct conformations of the single water complex H 2 0»S0 2 were 
found to be bound. These are shown in Fig. 12 together with some of the 
minimum energy geometry information. The total energies for the complexes 
and the molecular fragments as well as the binding energies are given in 
Table 19. 

For conformation I the geometry is almost exactly the same as it was 
with the earlier calculation with the larger basis set and the binding en- 
ergy is slightly greater by 1.5 mH. For conformation II the complex is 
found to be a little more compact and the binding energy is 2.6 mH greater 
than with the larger basis set. It should be noted that the gradient calcu- 
lation did not freeze the geometries of the fragments. Conformation II 
should be viewed as a nascent sulfurous acid molecule (H 2 S0 3 ) since a 

hydrogen transfer and subsequent redistribution of electrons would give 
o 

II 

HOSOH the presumed structure of sulfurous acid. 

Conformations III and IV represent single and double hydrogen-bonded 
structures, respectively. In conformation III the out-of-plane angle of 
rotation of the nonhydrogen-bonded OH is very soft energetically and no 
angle is favored energetically by more than a few tenths of a mi 1 1 iHartree. 
Conformation V was suggested by the results of a semiempirical calculation 
(Ref. 20), but tne minimum energy geometry result here indicates a 
preference for a cis arrangement rather than the trans arrangement of the 
semiempirical calculation. 
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Two-water complexes were studied using the energy gradient program with 
the 4-31G basis set. The conformations are shown in Fig. 13. Conformations 
VI and VII are mirror image doublings of conformations II and V, respec- 
tively. The interatomic distance did not change significantly in going to 
the two-water complexes. The energies are given in Table 19. In both cases 
the binding energies of the two-water complexes are a little less than twice 
the binding energies of the one-water complexes indicating relatively little 
interaction between the water molecules and almost independent addition of 
the second water molecule. 

To examine the effect of basis set size on the binding energy and ge- 
ometry two additional SCF energy gradient calculations were performed on 
conformation V. In the first calculation a double zeta basis set with 
McLean exponents was employed (Ref. 31). The binding energy decreased by 
3mH, a rough measure of the superposition error in the poorer basis set. In 
the second calculation a set of sulfur 3d orbitals (a = 0.6) was added to 
the double zeta basis set. The polarization functions reduced the S0 2 

bond angle, decreased the SO 2 dipole moment, and lowered the binding en- 
ergy of the complex by 2.4 mH. The results of these two calculations are 
given in Table 19. 

In order to examine the energetics of adding multiple waters to S0 2 
and still retain the double zeta plus sulfur d orbital basis set, SCF 
calculations were run using a generally contracted basis set at fixed 
geometries predicted by the earlier gradient results. The primitive basis 
set was a Huzinaga sulfur (12s/8p) oxygen (9s/5p), hydrogen (4s) set 
generally contracted to <4s/3p>, <3s/2p>, and <2s>, respectively, 

(Refs. 7,28). This set was augmented by a single set of 3d functions 
(a = 0.6) on the sulfur atom. This basis set was similar to the previous 
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double zeta plus d set and the binding energy of conformation V obtained 
with this new basis set differed by only 0.2 mH from the previous result. 

The multiple-water complexes studied are shown in Fig. 13. The total 
energies and binding energies are given in Table 19. Conformation I with 
this basis set produces a binding energy which is 72% of the binding energy 
derived from the earlier 4-31G calculation. This reduction in binding 
energy is virtually identical to that obtained for conformation V. The 
two-, three-, and four-water complexes differ by only a few mH from the sum 
of the binding energies of the separate bound pairs. This assumes that the 
energies of the conformations II and IV (which were not run with the new 
basis) can be approximated by taking 70% of the binding energy obtained with 
the 4-31G basis set. The binding energy of the five-water complex is, 
however, significantly less than the sum of the binding energies of the 
separate bound-pairs. 

The implications of the additivity of the binding energies supports the 
notion that the bonding interaction is nearly purely electrostatic. Only 
when the waters were close enough to interact with each other, as in the 
five-water complex, did the binding energy fall significantly below the sum 
of the binding energies of the separate bound pairs. An interesting 
corollary to the approximate additivity of the binding energies occurred in 
the three-water complex which was slightly more stable than predicted from 
the binding energy additivities perhaps due to a favorable alignment of the 
water dipoles. 

Although a Morakuma component analysis (Ref. 32) of the energy 

components of the H2O - SO2 interaction was not made, the results obtained 
here suggest that tne analysis would follow that for the NH 3 *S0 2 complex 
(Ref. 23) as discussed by Kollman (Ref. 33). The charge transfer and 
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polarization affects are small and the minimum energy geometry is determined 
by a balance of the electrostatic energy and the exchange repulsion energy. 
However, unlike the NH 3 .S0 2 complex, the minimum energy geometry for the 

H 2 0-SC)2 complex (conformation V) is close to that predicted solely from 
consideration of dipole-dipole interaction. 

The H-bonded complex (conformation III) has a binding energy (7.5 mH = 
4.7 Kcal) which is virtually the same as that obtained for the HF*S0 2 com- 
plex (5Kcal) (Ref. 24). The minimum energy geometry complex (conformation 
V) has a binding energy (17.4 mH = 10.9 Kcal) very close to that obtained 

for the NH3*S02 complex (9.3 Kcal) (Ref. 23). It should be noted that 
inclusion of d orbitals in the sulfur basis set reduced the binding energy 

(10.4 Kcal -► 9.3 Kcal) of NH3«S02, an effect that was observed in the 
H 2 0. S0 2 calculations. The CNDO/2 results of Castleman gave a much larger 
binding energy (145.9 mH = 98.3 Kcal) and a different minimum energy 
geometry (Ref. 20) . 

In the analysis of his molecular beam electric deflection results 
Castleman remarks on the existence of species with large dipole moments that 
he attributes to H 2 0.S0 2 complexes (Ref. 22). The dipole moments for some 
of the complexes were calculated in this study and are reported in Table 20. 
These values are all overestimated particularly those from calculations not 
using sulfur d orbitals. 

The geometry optimization routine also produced force constants and 
fundamental frequencies for vibrational modes. However, the results did not 
converge to sufficient accuracy to predict the sizes and signs of the fre- 

quency shifts of SO 2 and H 2 O vibrations upon complexation. Better conv- 
ergence would allow comparison with the experimental IR spectroscopic re- 
sults of Tso and Lee who claim to have evidence for the existence of the 
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H 2 0.S0 2 complex (Ref. 21). 

There are two implications from these results for acid rain chemistry. 
The most widely accepted mechanism for the tropospheric, gas phase produc- 


tion of H 2 SO 4 is as follows: 

1 

O 3 + hv ■"0( D) 




+ 

0 2 

( 1 ) 

0( D) + H 2 0 -2H0 



( 2 ) 

HO + S0 2 + M -H0S0 2 

+ 

M* 

(3) 

HOSO 2 + 0 2 — H0 2 

+ 

S0 3 

(4) 

S0 3 + H 2 0 -H 2 S0 4 



(5) 


This ab initio study shows that an H 2 OSO 2 complex is predicted to be 
stable and therefore an alternative to steps ( 2 ) and ( 3 ) above is 
suggested: 

H 2 0 + S0 2 — H 2 O.S0 2 (2') 

1 

0( D) + H 2 O.SO 2 — OH + HOSO 2 (3') 

where the complex could act as its own energy sink negating the need for a 
third body. A comparison of the activation energies for the alternatives is 
needed to decide between them. 

The second implication for acid rain chanistry comes from the results 
of the multiple-water complexes. There appears to be a shell of bonding 
sites for H 2 0 around a central S0 2 . This suggests the formation of a 
nucleation center with additional waters outside the initial sphere added in 
a typical H-bonded fashion to waters already in the complex. The energetics 
are favorable but a discussion of the position of equilibrium would have to 
include an evaluation of the entropy changes for the process. The formation 

of an S 02 ~containing droplet would allow heterogeneous chemistry to 
proceed. 
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In summary an energy gradient program was used to locate stable struc- 
tures for H 2 OSO 2 complexes and SCF calculations were carried out to de- 
termine the binding energies of complexes with multiple water molecules. A 
4-31G basis set was used for most potential energy searches. More accurate 
oasis sets including a generally contracted basis set with d orbitals on 
the sulfur atom were used for geometry and binding energy verification. For 
one water complexes five different stable geometries were located, with 
binding energies between 4 and 11 Kcal mol- 1 suggesting a binding shell for 
H 2 0 around S0 2 . Very little charge transfer between S0 2 and H 2 o was 
present. Addition of more than one H 2 0 was found to be energetically 
favorable and approximately additive up to four added waters. An 
alternative mechanism for the tropospheric, gas phase production of acid 
rain is suggested which does not require a third body. A process for the 
formation of an S0 2 -containing water droplet is advanced. 
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Table 17. Total energies (Hartrees) and binding energies (mi 1 1 iHartrees in 
parentheses) of conformations I and II and the separated 
fragments showing the effects of configuration interaction. 



SCF 

CI1 

CI2 

Separated fragments 

-622.9726 

-623.1224 

-623.1793 

Conformation I 

-622.9786 

(-6.0) 

-623.1339 

(-11.45) 

-623.1539 

(+25.4) 

Conformation II 

-622.9824 

(-9.8) 

-623.1457 

(-23.3) 

-623.1918 

(-12.5) 
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Table 18. Mul liken population analysis of conformations I and II and 

separated fragments showing effect of configuration interaction. 



Atom 

SCF 

Cl 2 

Separated fragments 

S 

+ 1.20 

+1.16 


Oi 

-0.60 

-0.58 


0 3 

-0.72 

-0.71 


H 

+0.36 

+0.355 

Conformation I 

S 

+1.23 



Oi 

-0.62 



03 

-0.74 



H 

+0.37 


charge transfer 

0.004 -►SO 2 


Conformation II 

S 

+2.35 

+1.20 


Oi 

-2.12 

-0.60 


O 2 

-0.75 

-0.61 


0 3 

-0.23 

-0.74 


Hi 

+0.37 

+0.37 


h 2 

+0.38 

+0.38 


charge transfer 0.52-^S0 2 0.001-^S0 2 
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Table 19. SCF total energies (Hartrees) and binding energies (mi 1 1 iHartrees 


H2 0 

in parentheses) of S 02 *(H 20 ) n 
Geometry Optimized 

-75.9086 -76.0110 

complexes. 

-76.0110 

Geometry fixed 
-76.0121 

S0 2 

-546.3733 

-546,9719 

-547.1540 

-547.1755 

I 

-622.2894(7.5) 

- 

- 

-623.1930(5.4) 

II 

-622.2943(12.4) 

- 

- 

- 

III 

-622.2894(7.5) 

- 

- 

- 

IV 

-622.2902(8.3) 

- 

- 

- 

V 

-622.2993(17.4) 

-622.9973(14.4} 

-623.1770(12.0) 

-623.1998(12.2) 

VI 

-698.2120(21.5) 

- 

- 

- 

VII 

-698.2236(33.1) 

- 

- 

-699.2212(21.5) 

VIII 

- 

- 

- 

-775.2437(31.9) 

IX 

- 

- 

- 

-851.2560(32.1) 

X 

- 

- 

- 

-927.2703(34.3) 

Basis 

Set A 

B 

C 

D 


A: 4.31G 

B: Double Zeta 

C: Double Zeta + Sulfur d 

D: S(12s/8p)/<4s/3p7, 0(9s/5p)/<3s/2p>, H(4s)/<2s> 
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Table 20. Dipole moments (in Debyes) for fragments and selected complexes 
from SCF geometry optimizations. 


H 2 0 

2.49 

2.53 

2.53 

S0 2 

- 

3.44 

2.67 

III 

6.10 

- 

- 

IV 

6.35 

- 

- 

V 

2.95 

2.92 

2.51 

Basis set 

A 

B 

C 
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h 2 

I II 


Figure 11. Conformations of S0 2 *H 2 0 complexes. 
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0 


0 



III IV 



V 

Figure 12. Conformations of SC^-I^O complexes including some minimum 
energy geometry results. 
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ADDENDUM 


Other projects under study during the grant period included some work 
on computational technique to study the above described systems. The most 
significant of these was the development and implementation of a method to 
calculate transition moments between non-orthogonal molecular orbitals by 
the use of corresponding orbitals. The details of the theory with some 
application have been published and a reprint of the article is attached 
with this report. A theoretical method was implemented to compute the 
vibrational frequencies, the overlaps between vibrational functions in a 
given electronic state and to determine the vibrational transition moments 
between various states. 

Other efforts included interfacing the BIGMOLI system of programs with 
SCF and Cl programs modified by other researchers at the Langley Research 
Center. Also the energy graduant integral program package received from the 
Quantum Chemistry Program Exchange was implemented at the Langley Research 
Center Computer Complex. 

Initial calculations on the triplet state potential energy surface 

H0 2 + HO + H 2 0 + 0 2 

were carried out using the energy gradient integral program package. The 
aim of these calculations was to study the stability and energetics of 
HO 2 OH complex. These calculations were postponed until the MCSCF procedure 
nas been incorporated into the energy gradient integral program package. 
Investigations of the lowest singlet state chemically bonded and hydrogen 
bonded intermediates in this reaction have been pursued in other efforts at 
the Langley Research Center. 
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The lower electronic states of CIOO: A computational investigation 
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Eight doublet and eight quartet states of CIOO were investigated by ab initio Cl techniques. The 
potential energy surfaces of the four lowest energy doublet states of both A ” and A ' symmetry 
indicate that only the 1 2 A * state is bound. In contrast to the model provided by the H0 2 radical, 
all of the excited doublet states investigated were repulsive with respect to dissociation to Cl + 0 2 
and metastable or bound with respect to dissociation to CIO + O. The transitions to the excited 
states investigated span the visible and near UV spectral regions, but the transition moments 
indicate that they are very weak. Since the photolysis products are the same as those of the rapid 
thermal dissociation, photolysis is not expected to be an important atmospheric process. The soft 
bending potential for the 1 2 A ' state and the shape of the 1 4 A " state in the entrance channel of the 
CIO + O — ► Cl + 0 2 reaction provide a qualitative explanation for the underprediction of the low 
temperature reaction rate by previous trajectory calculations. 


I. INTRODUCTION 

Destruction of stratospheric ozone by reaction ( 1 ), 

Cl + 0 3 ^C10 + 0 2 (1) 

is accompanied by the side reactions 

Cl + 0 2 + M X CIOO + M, (2) 

and 

CIO + CIO X CIOO + Cl. (3) 

The ratio of the rates of ( 1 ) and (2) and of the concentra- 
tion of ozone to that of molecular oxygen are such that reac- 
tions (1) and (2) occur at approximately equal rates in the 
upper atmosphere where ozone destruction by chlorine from 
anthropogenic sources is an important concern. Reaction (3) 
is not considered a major factor in the atmosphere, but may 
be important in laboratory investigations of the chemistry of 
ozone and chlorine. 1 In spite of its high rate of formation by 
reaction (2) the CIOO molecule is not considered very impor- 
tant in atmospheric chemistry because it rapidly dissociates 
by the reverse reaction. If the CIOO radical undergoes rapid 
reactions, however, these additional reactions could be im- 
portant to the atmospheric chemistry of chlorine and/or 
other trace species. A reaction involving the oxidation of 

al National Research Council Associate. Present address: Ballistics Re- 
search Laboratory, ORDAR-BLL, Aberdeen, MD 21005. 
b ‘ Present address: NASA-Ames Research Center, STS 230-3, Moffett 
Field, CA 94035. 


another atmospheric species, such as the oxidation of CO by 
CIOO proposed by Prinn 2 to explain the CO/C0 2 ratio on 
Venus, could play a role; but no rapid reaction of that type 
has been observed. Photolysis is another reaction which, de- 
pending upon the active wavelengths and products, could be 
important. Johnston, Morris, and Van den Bogaerde 3 ob- 
served a portion of the ultraviolet spectrum of CIOO by the 
use of modulation spectroscopy, but did not determine the 
products. In addition, Eachus et a/. 4(a| tentatively assigned 
an absorption feature observed at 380 nm in the photolysis 
products of C10 2 (in H 2 0/H 2 S0 4 matrices at 77 K and in 
crystalline KC10 4 ) to the CIOO radical. The production of 
0( ! i) ) is energetically possible at the wavelengths of John- 
ston et a/.’s observations and 0( 3 P) could be produced at 
longer wavelengths, which would not be subject to screening 
by ozone absorption in the atmosphere. Eachus et a/. 4(a| ob- 
served thermal isomerization of the CIOO radical to the 
symmetric C10 2 . Photoisomerization from CIOO to OCIO 
could also be significant, but only the reverse process was 
observed in Arkell and Schwager’s 4(b) matrix isolation inves- 
tigations of the infrared spectrum of CIOO. Other experi- 
mental investigations of CIOO have provided information 
on its heat of formation and rates of formation and dissocia- 
tion. 5 In a theoretical study, Gole and Hayes 6 evaluated the 
energy of CIOO relative to that of OCIO. They used the re- 
stricted SCF method for OCIO and CIOO - , but had to use 
the orbitals of the latter calculation for neutral CIOO be- 
cause of convergence difficulties. 

The present investigation is focused on providing 
further information on the photolysis spectrum and pro- 
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ducts: 

CIOO -f hv — ^products, (4) 

in order to provide for an assessment of the atmospheric 
importance of this process. Information on the lower states 
of CIOO is also of interest with regard to the reaction 

ClO + 0-*Cl + O 2 (5) 

which is an important element of the catalytic cycle for 
ozone destruction. A theoretical investigation 7 of this reac- 
tion relied on several assumptions regarding the lower states 
of CIOO and apparently underestimated the reaction rate. 
Recent experimental investigations 8 obtained approximate- 
ly the same rate at 298 K, but gave significantly different 
temperature dependencies. The results from the investiga- 
tion of the states of CIOO, which is an intermediate in Eq. (5), 
will be useful for understanding this reaction and the under- 
prediction of the rate by the trajectory calculations. 

Finally, understanding the electronic states of CIOO 
will provide a model in addition to that of H0 2 , for under- 
standing the range of photochemical behavior to be expected 
of peroxy compounds. Specifically, the states of FOO which 
may be of significance in atmospheric fluorine chemistry 
should be closely related to those of CIOO. The size of the 
CIOO molecule, the number of electronic states of interest, 
and the nature of its photochemical properties make this 
problem a significant computational task. The methods used 
to accomplish this task are discussed in the following sec- 
tion. The results obtained are presented and discussed in 
subsequent sections. 

II. COMPUTATIONAL METHODS 
A. Basis sets 

The (9s5/>) Cartesian Gausian basis of van Diujenveldt 9 
for oxygen and the {\2s9p) set of Huzinaga 10 for chlorine 
were generally contracted 1 1 to {3s2p) and (4s 3 /?), respectively, 
and augmented by a set of s and p bond centered functions 12 
to provide for polarization. The exponents of the bond func- 
tion were 0.7 ( 5 ) and 0.6 ( p). 

The molecular orbitals were taken from a limited “pairs 
only” MCSCF treatment of theX 2 A " ground state in which 
the SCF configuration (CSF) 

^ SCF = . . .(10a') 2 (11a') 2 (12a') 2 (13a') 2 {2a"') 2 (3a") 2 4a" 

was augmented by CSF’s in which the ct( 0-0) and <j(C 1-0) 
bonding orbitals were replaced by their respective antibond- 
ing orbitals. Examination of Cl treatments utilizing orbitals 
from more extensive pairs only MCSCF treatments and 
from MCSCF treatments of the 1 2 A ' state indicated that the 
above choice is superior. In the SCF description of theX 2 A " 
state, the 10a' orbital is a lone pair on the terminal oxygen, 
11a' is the 0-0 bonding orbital, 12a' is the Cl lone pair 
orbital and 13a' is the Cl-O bonding orbital. The out-of- 
plane orbitals 2A ", 3a", and 4a" are somewhat localized on 
Cl, the central O atom, and the terminal O atom, respective- 
ly* 

Preliminary configuration interaction (Cl) calculations 
indicated that all of the configurations of importance arise 
from rearrangements of the electrons within the occupied 


orbitals and excitations to the 14a' and 15a' orbitals which 
are the Cl-O and 0-0 antibonding orbitals, respectively. 
The results of these calculations further indicated that exci- 
tations from the molecular orbitals arising from core orbitals 
and psuedo core orbitals arising primarily from the Cl (3/*) 
and O (25 ) components were not important to any of the 
states of interest. Although the most important CSF’s for 
describing the lower states of either symmetry were the same 
at widely varying points on the potential energy surfaces, a 
relatively larger number of CSF’s were of minor importance 
at some points on the surface. At a given point on the poten- 
tial energy surfaces between 20 and 50 CSF’s made at least a 
1% contribution to one of the four lowest doublet states. 
Supplementary calculations on the effect of varying the ref- 
erence list across the surface indicated that only minor 
“noise” would result. This is in agreement with recently pub- 
lished investigations. 13 The reference CSF’s actually used at 
each geometry were selected from “frozen core” full valence 
Cl calculations on the lowest four roots of each symmetry. 

A list of target CSF’s was generated by including all 
single and double excitations from the reference list subject 
to the restrictions (1) that the true core orbitals remain fully 
occupied, (2) that only one hole could exist in the orbital 
space arising from the chlorine 3 P and oxygen 25 atomic 
orbitals, and (3) that only one electron could occupy the ex- 
ternal (nonvalence) portion of the orbital space. 

A final CSF list was prepared by selecting CSF’s from 
the target list according to their A k contributions 1 3(b} to the 
first four roots using a cumulative threshold of 60 mhartree. 
In most calculations all single excitations with respect to the 
most important CSF for each of the four states of interest 
were retained without regard for their energy contributions. 

The perturbation theory estimates of the energy contri- 
butions, to each root, for all rejected CSF’s were added to the 
raw eigenvalues in order to obtain comparable energies for 
the different states and for the same state at different geome- 
tries. Extensive calculations were carried out using smaller 
thresholds to determine the extent of the potential energy 
surface noise introduced by the selection and energy adjust- 
ment procedures. The effect of selection threshold on the 
convergence of dipole moments and transition moments was 
also investigated. The potential energy noise was found to be 
about 10% of the cumulative threshold in amplitude. This 
agreed with earlier unpublished results. 14 Although the 
computed dipole moments and transition moments were 
quite small, the results obtained using the procedures dis- 
cussed above led to values which were converged to within 
40% of the values obtained using small thresholds. The over- 
all results of the calculations which were carried out to assess 
the accuracy of the present methods was to insure that the 
results would provide sufficient accuracy for the major goals 
of this study. In addition to the above Cl calculations, a few 
additional calculations were carried out on theX 2 A " ground 
state using a smaller (5 mhartree) threshold and natural orbi- 
tals obtained from the above calculations. The reference 
wave function for the natural orbital calculations included 
all CSF’s which contributed 0.25% or more to a (frozen core) 
full valence Cl wave function. The target CSF list for the 
natural orbital Cl was formed using the prescription de- 
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TABLE I. Energies of 2 A 0 and 2 A ' states* 


Calc# 

R (CIO) 

*(00) 

ClOO 

1 W 

mm 

3 2 A* 

4 2 A ’ 

1 2 A’ 

EH 

3 2 A ' 

4 2 A ' 

1 

3.6 

2.0 

113.6 

- 0.2019 

- 0.0828 



-0.1251 

- 0.0867 

- 0.0498 

- 0.0393 

2 

3.6 

2-2 

113.6 

-0.2711 

-0.1585 

-0.1379 

-0.1285 

-0.1917 

-0.1706 

-0.1244 

-0.1242 

3 

3.6 

2.48 

113.6 

-0.2881 

-0.1808 

-0.1693 

-0.1489 

-0.2211 

- 0.2077 

-0.1563 

-0.1489 

4 

3.6 

2.6 

113.6 

- 0.2809 

-0.1762 

-0.1668 

-0,1494 





5 

3.6 

2.8 

113.6 

- 0.2636 

-0.1689 

-0.1420 


- 0.2247 

-0.1833 

-0.1344 

-0.1352 

6 

3.6 

3.1 

113.6 

- 0.2378 

-0.1600 

HHI 

- 0.0966 

-0.2157 

-0.1628 

-0.1237 

-0.1078 

7 

3.6 

3.5 

113.6 


-0.1415 

BH9 

-0.1299 

-0.1927 

-0.1542 

-0.1324 

-0.0801 

8 b 

3.6 

4.2 

11 3.6" 

-0.1713 

-0.1538 

-0.1509 


-0.1817 

-0.1621 

-0.1452 

-0.1157 

9 b 

3.6 

5.0 

1 1 3.6 b 


-0.1596 

-0.1566 

- 0.0802 

-0.1631 

-0.1617 

-0.1593 

- 0.0789 

10 

3.1 

2.48 

113.6 

- 0.2717 

-0.1134 

-0.0648 

- 0.0254 





11 

3.3 

2.48 

113.6 


-0.1485 

-0.1185 

- 0.0948 

-0.2192 

-0.1657 

-0.1134 


12 

3.45 

2.48 

113.6 

- 0.2848 

-0.1642 

-0.1485 

-0.1462 





13 

3.75 

2.48 

113.6 

— 0.2881 

-0.1969 

-0.1815 

-0.1776 

-0.2219 

-0.1973 

-0.1684 

-0.1657 

14 

3.90 

2.48 

113.6 

- 0.2878 

-0.2124 

- 0.2048 


- 0.2303 

-0.2160 

-0.1939 


15 

4.3 

2.48 

113.6 

— 0.2910 

- 0.2423 

- 0.2355 

- 0.2902 





I6 b 

5.0 

2.48 

113.6 b 

- 0.2862 

-0.2718 

-0.2458 

- 0.2358 

- 0.2470 

- 0.2450 

- 0.2334 

-0.2115 

I7 b 

7.0 

2.48 

113.6 

- 0.2829 

- 0.2794 

- 0.2449 

- 0.2338 





18” 

25.0 

2.40 

180.0 

-0.2816 

- 0.2809 

- 0.2408 

- 0.2383 





19 

3.75 

2.48 

80.0 

- 0.2206 

-0.1943 

-0.1732 

-0.1214 

-0.2019 

-0.1639 

-0.1342 

-0.1255 

20 

3.75 

2.48 

95.0 

- 0.2704 


-0.1792 

-0.1560 

- 0.2232 


-0.1642 

-0.1478 

21 

3.75 

2.48 

105.0 

- 0.2847 

-0.1989 

-0.1808 

-0.1736 

- 0.2277 

- 0.2099 


-0.1602 

22 

3.75 

2.6 

113.6 


-0.1914 

-0.1805 

-0.1707 





23 

3.75 

2.48 

122.0 

- 0.2863 

-0.1977 

-0.1782 

-0.1763 

- 0.2259 

-0.2111 

-0.1722 

-0.1712 

24 

3.75 

2.48 

140.0 

- 0.2722 

- 0.2010 

-0.1835 

-0.1692 

-0.2188 

- 0.2028 

-0.1711 

-0.1646 

25 

3.75 

2.48 

160.0 

- 0.2493 

- 0.2076 

-0.1880 

-0.1593 

-0.2131 

-0.1919 

-0.1730 

-0.1506 

26 

3.75 

2.48 

180.0 

- 0.2221 

-0.1912 

-0.1711 

-0.1090 

- 0.2070 

-0.1961 

-0.1807 

-0.1562 

27 

3.75 

2.3 

180.0 

- 0.2335 

-0.2016 

-0.1961 

-0.1458 





28 

3.75 

2.8 

180.0 

— 0.2018 

-0.1858 

-0.1738 

-0.1381 





29 

3.75 

2.1 

160.0 

-0.2110 

-0.1616 

-0.1445 






30 

3.75 

2.3 

160.0 

-0.2464 

-0.2021 

-0.1849 

-0.1464 





31 

3.75 

2.8 

160.0 

-0.2126 

-0.1742 

-0.1454 

-0.1364 





32 

3.75 

3.1 

160.0 

-0.1827 

-0.1494 

-0.1421 

-0.1177 





33 

3.75 

2.36 

113.6 

-0.2906 

-0.1968 

-0.1819 

-0.1721 





34 

3.2 

3.5 

113.6 

-0.2194 

-0.1560 

-0.1477 

-0.1362 





35 

3.4 

3.5 

113.6 

-0.2161 

-0.1497 

-0.1438 

-0.1267 





36 

3.3 

2.3 

113.6 

-0.2710 

-0.1237 


-0.0817 





37 

3.3 

2.6 

113.6 

-0.2781 

-0.1282 

- 0.0957 

- 0.0683 





38 

3.6 

2.48 

105.0 

- 0.2832 


-0.1663 

-0.1504 





39 

4.1 

2.48 

113.6 

- 0.2891 

- 0.2279 

- 0.2260 






40 

4.5 

2.48 

113.6 

- 0.2896 

- 0.2549 

- 0.2399 

-0.2190 






a Energies in hartrees relative to — 609.0 hartree. Although energies are given to four decimal places the uncertainty introduced by the selection and 
extrapolation procedure occurs in the third decimal place (see the text). 
b Forced retention of single excitation was not used in these calculations. 



FIG. 1 . Cl energies of' 2 A " states as a function of R (CIO). The ClOO angle is 
constant at 1 13.6 deg and R (OO) equals 2.48 bohr. 



R(O-O), BOHR 

FIG. 2. Cl energies of the 2 A " states as a function o(R (OO). The ClOO angle 
is 1 13.6 deg and R (CIO) is 3.6 Bohr. 
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ORIGINAL PACS - 

OF POOR QUALITY 

scribed above. 

III. RESULTS 

The energies obtained for the eight doublet states treat- 
ed are presented in Table I for each of the geometries investi- 
gated. The natures of the potential energy surfaces are dis- 
cussed below. Some of the data from Table I is duplicated in 
figures for the sake of clarity. 

A. The 2 A * states 

Slices through the potential energy surfaces of the four 
lowest 2 A " states along the CIO and OO bond stretching 
coordinates are shown in Figs. I and 2. 

1. The X 2 A" state 

The potential energy surface along the CIO bond 
stretching coordinate (Fig. 1) is broad and shallow as would 
be expected for a weak bond. The surface exhibits a weak 
double minimum character but this is probably due to CSF 
selection noise (vida supra). 

The electronic binding energy relative to that of Cl + 0 2 
products is calculated to be approximately 0.01 hartree 
(6 kcal/mol) in agreement with the experimental value of 
7.6 kcal/mol. However, the present treatment is focused on 
providing a description of several excited states and this re- 
sults in a rather large uncertainty ( ± 100%) in the comput- 
ed D e . The natural orbital Cl, which should be somewhat 
more accurate, gave similar results (Table II). We, therefore, 
estimate that D e = 7.5 ± 4 kcal/mol. 

The shape of the potential along the 0-0 stretching co- 
ordinate (Fig. 2) appears much more conventional and the 
curvature near the bottom of the well agrees with Arkeil and 
Schwager’s experimental force constant for the 0-0 stretch. 
The calculations indicate that ground state is bound by 3.46 


TABLE II. Cl energies of the X 2 A * state using approximate natural orbi- 
tals. 3 


R{ CIO) 

R{ OO) 

<cioo 

E 

3.1 

2.48 

113.6 

~ 609.2673 

3.3 

2.48 

113.6 

- 609.2796 

3.45 

2.48 

113.6 

- 609.2841 

3.6 

2.48 

113.6 

- 609.2867 

3.75 

2.48 

113.6 

- 609.2874 

3.9 

2.48 

113.6 

- 609.2844 

4.1 

2.48 

113.6 

- 609.2860 

4.3 

2.48 

113.6 

- 609.2848 

4.5 

2.48 

113.6 

- 609.2829 

5.0 

2.48 

113.6 

- 609.2783 

7.0 

2.48 

113.6 

- 609.2742 

25.0 

2.40 

180. 

- 609.2774 

3.60 

2.20 

113.6 

-609.2681 

3.60 

2.48 

113.6 

- 609.2867 

3.60 

2.60 

113.6 

- 609.2745 

3.75 

2.20 

105. 

- 609.2755 

3.75 

2.48 

105. 

- 609.2847 

3.75 

2.60 

105. 

- 609.2785 

3.75 

2.36 

113.6 

- 609.2887 

3.75 

2.48 

122. 

- 609.2892 

3.60 

2.60 

105. 

- 609.2831 


a Bond lengths in bohr, angles in degrees, and energies in hartrees. 


eV relative to dissociation to CIO and O. This value shoul 
be reduced by — 0.8 eV to account for relaxation of the Cl< 
bond length at dissociation. 15 The resulting value, whe 
compared to the difference between the dissociation energk 
of CIO and 0 2 gives — 7 kcal/mol for the binding energy c 
CIO relative to CIO + 0 2 in (fortuitously) good agreemer 
with the value obtained by direct calculation. The groun 
state dipole moment is relatively small ( — 0.1 a.u.) in th 
equilibrium region (Table III). At shorter bond lengths th 
charge transfer and resulting dipole moment increases. 

2. The 2 A * excited states 

The most important feature of the three 2 A " excite* 
states treated is that they are strongly repulsive with respec 
to dissociation into Cl and 0 2 fragments (Fig. 1). All three o 
the 2 A " excited states lie in a narrow band in the equilibriun 
region. The vertical excitation energies are calculated to fal 
in the range 2.0 to 3.8 eV with an uncertainty of —0.3 eV 
The first excited state dissociates to the ground state asym 
pote, ( 2 P) Cl -f ~ )0 2 , while the second and third excitec 
states apparently yield ( l zl g )0 2 . 

The behavior of the 2 A " states as a function of the O 
bond length is shown in Fig. 2 for fixed values of R (CIO) anc 
< ClOO. A barrier to dissociation of the 2 2 A " state tc 
CIO 4- O is predicted to occur at — R (OO) = 3.5 bohr. If the 
R (Cl-O) bond length is relaxed at R (0-0) = 3.5, the barrier 
is predicted to be 0.7 eV. In addition, if the CIO potential 
energy of Arnold, Whiting, and Langhoff 15 is used to esti- 
mate the energy lowering which would be obtained by relax- 
ing the CIO bond length at infinite R (0-0), the 2 2 A " state is 
found to be slightly metastable. Due to the uncertainties in 
the present energies, this prediction is not definitive. The 
second and third excited 2 A " states also exhibit a potential 
well with respect to R (0-0) in the equilibrium region. 

The third excited 2 A " state has a maximum at 
~R (0-0) = 3.1 bohr which, on the basis of the transition 
moment to the ground state (vide infra), is apparently due to 
an avoided crossing with a higher state. An avoided crossing 
occurs between the 3 2 A " and 2 2 A " states near R (O- 


TABLE III. Dipole moment of the X 2 A ” state. 3 


R (CIO) 

(OO) 

ClOO 

M, MSCSF 

M, Cl 

3.75 

2.48 

95 

0.045 

0.124 

3.75 

2.48 

105 

0.046 

0.111 

3.75 

2.48 

113.6 

0.051 

0.105 

3.75 

2.48 

122 

0.056 

0.110 

3.6 

2.0 

113.6 

0.358 

0.579 

3.6 

2.2 

113.6 

0.071 

0.279 

3.6 

2.48 

113.6 

0.091 

0.027 

3.6 

2.6 

113.6 

0.120 

0.060 

3.6 

2.8 

113.6 

0.145 

0.147 

3.6 

3.1 

113.6 


0.227 

3.6 

3.5 

113.6 

0.151 

0.197 

3.3 

2.48 

113.6 

0.202 

0.327 

3.45 

2.48 

113.6 

0.145 

0.109 

3.9 

2.48 

113.6 

0.037 

0.147 

5.0 b 

2.48 

113.6 

0.041 

0.107 


3 Dipole moment in atomic units, bond lengths, in bohr, and angles in de- 
grees. 

b Single excitations which were rejected during configuration selection were 
not retained in this calculation. 
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R(Ci-O), BOHR 

FIG. 3. Cl energies of the 2 A ' states as a function of R (CIO). The CIOO angle 
and R (OO) bond length are the same as those given in Fig. 2. 


O) = 3.5 bohr. The character of the 2 A " excited states 
changes drastically for R (0-0)>3.5 bohr. At R (O-O) = 5 
bohr the first three 2 A " states are converging to ground state 
products while the energy of the 4 2 A " state indicates that it 
dissociates to yield 0( l D ). 

The bending potentials for the 2 A " excited states 
(Table I) exhibit an orbital basis artifact for the excited states. 
However, aside from this artifact, the two lowest 2 A " excited 
state bending potentials are relatively flat. Several calcula- 
tions (Table I) were carried out to explore the possibility of a 
low energy channel for dissociation to CIO + O. No lower 
energy channel for dissociation of the excited states to these 
products was revealed by these calculations. 

B. 2 A ' states 

Slices through the 2 A ' potential surfaces along the 
R (CIO) and R (OO) coordinates are shown in Figs. 3 and 4. 
As with the 2 A " excited states, all four of the 2 A ' states were 
repulsive with respect to dissociation into Cl and 0 2 , had 
bound regions with respect to dissociation into CIO + O, 



FIG. 4. Cl energies of the 2 A ' states as a function of R (OO). The CIOO angle 
and CIOO bond length are the same as those given in Fig. 3. 



2.0 3.0 4.0 5.0 

R(O-O), BOHR 


FIG. 5. Cl energies of 4 A ' states as a function of R (OO). The CIOO angle 
and CIO bond lengths are the same as those given in Fig. 3. 

and exhibited an artifact “bump 1 ’ in the bending potentials 
due to the orbital bias. This bump is evident in the CIO 
stretch potential for the 2 2 A ' state (Fig. 3) at R (CIO) — 3.75 
bohr. Neglecting the artifact, the bending potentials (Table I) 
of the 2 A ' states appear relatively flat. The behavior of the 
potentials at small angles may be indicative of a low energy 
path between the excited states of CIOO and OCIO. The first 
two 2 A ' states lie below the 2 2 A " state. The existence of one 
low lying 2 A ' state is expected on the basis of the near in- 
frared spectra of H0 2 and other peroxy compounds. The 
first three 2 A ' states dissociate to ground state products on 
the R (CIO) coordinate and the 4th 2 A ' state produces 
( l A g )0 2 . On the R (OO) coordinate the first three states yield 
ground state products while the fourth dissociates to ground 
state CIO and ( l D )0. However, the third and fourth states 
are predicted to have intersections on both the inner and 
outer limbs so the fourth state may be effectively predisso- 
ciated to ground state products. 

In summary, the 2 A ' states investigated are expected to 
lead to absorption in the 1.5 to 5 eV region. The quantum 
yield should be unity for the entire spectral range. The pho- 



FIG. 6. Cl energies of the 4 A ' states as a function of R (OO). The CIOO angle 
and the CIO bond length are the same as those given in Fig. 3. 
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TABLE IV. Transition moments between X 2 A * and the doublet excited states,* 


R (CIO) 

Geometry 

(OO) 

CIOO 

2 2 A ' 

Upper state 
3 2 A” 

4 2 A " 

1 2 A m 

2 2 A ' 

3 y A ' 

4 2 A ' 

3.1 

2.48 

113.6 

0.433 

0.116 

1.035 





3.3 

2.48 

113.6 


0.112 

0.066 

0.090 

0.148 

0.030 

0.308 

3.6 

2.48 

113.6 

0.063 

0.081 

0.035 

0.014 

0.002 

0.007 

0.069 

3.75 

2.48 

113.6 

0.062 

0.053 

0.025 

0.012 

0.010 

0.046 

0.030 

3.90 

2.48 

113.6 

0.042 

0.032 

0.023 

0.005 

0.003 

0.040 

0.002 

5.0 

2.48 

113.6 

0.032 

0.032 

0.010 

0.044 

0.019 

0.007 


3.6 

2.0 

113.6 

0.045 

0.050 

0.058 

0.005 

0.005 

0.074 

0.005 

3.6 

2.7 

113.6 

0.037 

0.037 

0.043 





3.6 

2.8 

113.6 

0.080 

0.078 

0.008 

0.006 

0.0 11 

0.047 

0.036 

3.6 

3.1 

113.6 

0.073 

0.050 

0.276 

0.005 

0.013 

0.030 

0.021 

3.6 

3.5 

113.6 

0.010 

0.038 

0.041 

0.003 

0.008 

0.001 

0.025 

3.75 

2.48 

95 

0.022 

0.255 

0.139 

0.014 

0.003 

0.043 

0.019 

3.75 

2.48 

105 

0.042 

0.123 

0.063 

0.0 11 

0.003 

0.050 

0.017 

3.75 

2.48 

122 

0.060 

0.036 

0.021 

0.007 

0.001 

0.024 

0.002 


* Transition moments in atomic units. 


tolysis products are predicted to be entirely Cl -F 0 2 in the 
visible and near UV region and predominately these pro- 
ducts in the higher energy region ( > 3.5 eV) where the chan- 
nel for CIO + O production is open. 


C. Quartet states 

The behavior of the four lowest^ " and A ' quartet states 
were investigated along the R (OO) coordinate. The results, 
which are given in Figs. 5 and 6, show the first two 4 A " states 
and the lowest energy 4 A ' state to be bound with respect to 
dissociation into CIO -F O. The remaining five quartet states 
are repulsive with respect to dissociation into CIO -F O. The 
three lowest states of each symmetry dissociate to ground 
state products while the fourth state in each symmetry ap- 
parently dissociates to higher states of the products. The 
repulsive quartet states predissociate the fourth 2 A " and 2 A ' 
states. Since the latter states are repulsive with respect to 
dissociation to Cl 0 2 , the quartet states are not expected 
to be important to the photodissociation of any of the states 
investigated here, but many may be significant for O -f CIO 
collisions. In particular, reactive scattering on the 1 4 A " sur- 
face may be responsible, in part, for the difference between 
the theoretical 7 and experimental 8 results for the rate of the 
reaction (vide infra). 

D. Transition moments 

The transition moments for the 2 A " and 2 A ' excited 
states are given in Table IV for several geometries. The tran- 
sition moments are quite small near the ground state equilib- 
rium geometry and become even smaller at larger values of 
R (CIO) as would be expected for excited states arising from 
the ground state of Cl and the ground and V A g states of 0 2 . 
The calculated values of R (CIO) — 5.0 bohr are a little larger 
than they should be because forced retention of single excita- 
tion was not used at that geometry. 

The transition moments become somewhat larger for 
CIO bond lengths and bond angles smaller than the values 
for the ground state equilibrium. The larger value obtained 
for the transition moment to the 4 2 A "state of R (O-O) = 3.1 


bohr is indicative of mixing between the state and a higher 
state of the same symmetry and may explain the behavior of 
the 4 2 A " energy in that region. 

These results indicate that the absorption of CIOO will 
be quite weak in the visible region. The intensity of the solar 
spectrum in the visible region will substantially enhance the 
importance of this weak absorption with regard to atmo- 
spheric processes. However, since the products of photodis- 
sociation at longer wavelengths are Cl and 0 2 , the net effect 
of photodissociation will be to slightly enhance the total rate 
of dissociation to these products; a process which is dominat- 
ed by the high rate of thermal dissociation. 

IV. DISCUSSION 
A. Photochemistry of CIOO 

The repulsive nature of the lower doublet states of CIOO 
on the CIO coordinate should dominate the photochemistry 
for wavelengths in the visible and near UV. The absorption 
strength at shorter wavelengths will involve states that were 
not included in the present investigation. The calculated ver- 
tical excitation energies for the states included in the present 
treatment range from approximately 1.8 to 3.2 eV. The ex- 
tent of the energy range of the absorption system due to these 
states will depend on the spatial extent of the nonzero ampli- 
tude of the lower vibrational levels of the ground state. Al- 
though the present ground state results are not sufficiently 
accurate to establish the latter, the system probably extends 
from between 4 and 5 eV down to about 1.5 eV. Although it 
is possible to estimate the spectrum on the basis of harmonic 
oscillator vibrational functions chosen to reproduce the ex- 
perimental CIO stretching frequency, such an estimate could 
err significantly near the absorption edges, particularly if the 
actual potential is quite anharmonic. Since the products of 
photolysis are the same as those of the rapid thermal disso- 
ciation, photolysis of CIOO should not have a significant 
effect in atmospheric chemistry. 

Arkell and Schwager 4,at observed photo-isomerization 
of OCIO to CIOO. The absence of significant repulsion in the 
small angle portion of the potential energy surfaces of several 
excited states appears to provide a mechanism for the photo- 
isomerization. The repulsive nature of the CIOO excited 
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states explains these authors failure to observe the reverse 
process. Arkell and Schwager suggested that the loss of 
ClOO observed at temperatures (above 20 K) high enough 
for the Cl atoms to be mobile might have been due to photo- 
lysis (to Cl 4- 0 2 ) by the infrared source in their spectropho- 
tometer. At lower temperatures the lattice cage effects 
would lead to immediate recombination of the photolysis 
products. 

The remaining observation by Arkell and Schwager 
which is of interest here is their assignment of infrared lines 
at 435 and 1414cm~‘ toathird isomer of C10 2 . Their results 
indicated that this species was produced by the photolysis of 
OCIO and converted to ClOO upon warming above 4 K. 
They assigned this species as a geometrical isomer of ClOO 
because the isotopic splittings indicated that it was unsym- 
metrical and its lifetime at 20 K was too long to support its 
assignment to a metastable electronic state. The present re- 
sults suggest two possibilities for this species: (1) that the 
double minimum nature calculated for the ground state is 
real and that the metastable isomer corresponds to the high- 
er energy of the two minima; and (2) that there exists a metas- 
table region in the small ClOO angle portion of the 1 2 A ' 
surface. As noted earlier, the shallow double minimum in 
the ground state is probably an artifact of the computational 
techniques. In order for the latter possibility to satisfy all of 
the experimental constraints, the metastable l 2 A ' state ge- 
ometry would have to correspond to a distorted ring state 
and would have to lie below the 1 2 A " state in that conforma- 
tion. Although the present calculations did not provide de- 
finitive evidence that the latter possibility is fulfilled, the 
shapes of the computed surfaces indicate that it might exist. 

We conclude this section by again pointing out that al- 
though the photochemistry of ClOO involves a large number 
of relatively low lying states, the repulsive nature common to 
the lowest doublet excited states is expected to be the domi- 
nant factor. 

B. Comparison of the states of ClOO and HOO 

It is useful to compare the results for ClOO, about which 
little was previously known, to those of H0 2 , H0 2 is the 
simplest peroxy compound and has been the subject of sever- 
al experimental 16 and theoretical 17 investigations. The simi- 
larities in the gross electronic structure of the 0 2 fragment 
and in the observed 18 low energy 1 2 A ' transition of 

several peroxy compounds has led to the states of H0 2 being 
used as a model for understanding the electronic spectra of 
other peroxy compounds. Some similarities also exist 
between the observed UV spectra of H0 2 and ClOO. 3 

Although the transition moments are larger at the con- 
formations which would lead to absorption in the quartz UV 
region investigated by Johnston et a/., 3 the Franck-Condon 
factors for these are small. Poorly converged results for high- 
er states (above the 4 2 A " and 4 2 A ' states) indicate that these 
higher states of ClOO are probably primarily responsible for 
the absorption observed at —2500 A. The transition mo- 
ment obtained for the transition to the 4 2 A " state at the 
avoided crossing with the 5 2 A " at R (OO) = 3.1 bohr indi- 
cates that transitions to these higher states carry a larger 
intensity. 


The ClOO molecule also has a relatively low lying 1 2 A ' 
state, but the present results indicate that it is repulsive with 
respect to dissociation into Cl 4- 0 2 and it has not been ob- 
served experimentally. Also, the gross similarity of the UV 
spectra of ClOO and HOO is fortuitous. The UV absorption 
in HOO is due to transitions to the 2 2 A" and 2 2 A ' states 
which are metastable with respect to breaking the X-O bond 
and repulsive with respect to the 0-0 bond coordinate. 
These states lie at lower energies in ClOO and exhibit the 
opposite behavior with respect to the shape of the surfaces. 
The X 2 A " -+2 2 A " transition carries considerable intensity 
in HOO, but is quite weak in ClOO. The differences between 
the electronic states of these structurally similar molecules 
arise from several factors, but the major factor is the weak- 
ness of the CIO bond. This leads to a small splitting between 
the bonding and antibonding orbitals for this bond, lowering 
the energy of <7<C10) — ► a (CIO) excitation and the excita- 
tions from the lone pair orbital on the terminal oxygen to the 
CIO antibonding orbital. The low energy of the latter results 
in this transition being effective in back donating electron 
density in the a space to compensate for charge transfer in 
the pi orbital space. It is the latter which is responsible for the 
ionicity 17(cU7(d) of the 2 2 A " state of H0 2 . The magnitude of 
the transition moment between the lower 2 A " states in HOO 
is due to this ionicity, 17(d) which is largely absent in ClOO. 
The dipole moment of the 2 2 A " state of ClOO at tne ground 
state equilibrium geometry is only 0.98 D, whereas the di- 
pole moment of the 2 2 A " state of HOO is 4 D. 1 7(d) The tran- 
sitions from the lone pair electrons on Cl also contribute to 
the dense spectrum of states in the visible through quartz UV 
regions, compared to that of HOO. 

C. Other peroxy molecules 

The infrared spectrum 19 of FOO indicates that the FO 
bond is weak. Since it has the same number of valence elec- 
trons as ClOO and is structurally similar, its spectrum and 
photolysis products would be expected to be similar to those 
predicted here for ClOO. A long wavelength band 
(300</l<600 nm with — 420 nm) has been ob- 

served 20 ” 22 for FOO and the photolysis products (above 
—490 nm in argon matrix) are reported 22 to be F 4- 0 2 . A 
shorter wavelength band A max —205 nm has also been ob- 
served, but the photolysis products are now known. The 
theoretical results obtained here for the lower energy transi- 
tions in ClOO provide support for the experimental assign- 
ments of the observed long wavelength band to FOO and for 
assignment of the products in that band to F 4- 0 2 . 

The CO bond strength in CH 3 0 2 is apparently subject to 
a large uncertainty. If the 35 electrons of the Cl atom are 
included in the valence shell, CH 3 0 2 is formally iso-valence- 
electronic to ClOO, but there are significant differences 
between the lone pairs of Cl and the bonding orbitals of CH 3 . 
Experimental evidence 18 indicates that the low lying 1 2 A ' 
state (assuming C s symmetry) is bound in CH 3 0 2 and a short 
wavelength band has been observed. The former suggests 
that at least two states of CH 3 0 2 are similar to those of K0 2 , 
but the question regarding the photolysis products and of the 
existence of longer wavelength absorption are not resolved 
experimentally. These questions are significant for the atmo- 
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spheric chemistry of CH 3 0 2 in the relatively clean tropo- 
sphere. Photolysis may compete with reactions of CH 3 0 2 
with other trace species if it possesses absorption beyond 300 
nm and would be important if it results in products other 
than CH 3 4 0 2 . 

D. The CIO + O — ► Cl + 0 2 reaction 

This reaction, which closes the loop in the catalytic de- 
struction of ozone by Cl, was the subject of several experi- 
mental investigations 8 and a theoretical investigation by 
Jaffe. 7 The experimental investigations 8(a),8(b) gave opposite 
temperature effects and the theoretical rate was lower than 
the experimental values. 

Jaffe’s trajectory calculations on the CIO 4 O 
— * Cl 4- 0 2 reaction utilized several assumptions regarding 
the states of the CIOO intermediate. These included assum- 
ing that; (a) the topography of the ground state and the 1 2 A ' 
state are similar so that the ground state potential could be 
used for both states; (b) the lowest quartet states are unreac- 
tive; and (c) that the shape of the surfaces are such that the 
OCIO region is not accessible at low temperature. The pres- 
ent results indicate that the bending potential for the 1 2 A ' is 
quite flat so that assumption (a) would result in some under- 
prediction of the reaction rate. However, Jaffe discounted (a) 
somewhat by using effective bending force constants of 1/2 
and 1/4 the experimental 2 A " value for some of his calcula- 
tions. The flat bending potential for the 1 2 A ' state also sug- 
gests that assumption (c) may be inappropriate. 

The present results are not sufficiently complete to en- 
able more than speculation on the latter point. The observa- 
tions 4 of thermal isomerization and photoisomerization in 
certain matrices are also suggestive, but not definitive in this 
regard. The net effect of introducing a low energy channel to 
the OCIO surface is not clear since it would result in pro- 
cesses which would reduce the calculated rate (especially at 
high temperatures) 

O 4 - OC1 — CIOO — OCIO — O 4- CIO, 
as well as those which would increase the calculated rate 

O 4- CIO — ► OCIO — CIOO Cl + 0 2 . 

The present calculations indicate that the 1 4 A " state 
does not have a large barrier in the entrance channel so that 
assumption (b) should contribute to an underprediction of 
the rate by the trajectory calculations. Neither the exit chan- 
nel nor the bending portions of the surface were investigated 
* here, but it is relatively safe to assume that the exit channel 
has the same repulsive character as the doublet excited 
states. 

The present results appear to qualitatively explain the 
underprediction of the low temperature rate by Jaffe’s trajec- 
tory calculations. Because of the three surfaces contributing 
to the low temperature rate, the additional surfaces which 
would become important at high temperatures, and the pos- 
sibility that both OCIO and CIOO intermediates may be in- 
volved; the mechanism cannot be regarded as being simple. 
Some additional calculations on the potential energy sur- 
faces and more elaborate trajectory calculations would be 
required to resolve the uncertainties 22 regarding this reac- 
tion. 


E. Other reactions 

Although reactions between CIOO and other species 
were not investigated directly in the present work, the infor- 
mation on the bonding and electronic properties does facili- 
tate speculation on CIOO reactions. 

The weak CIO bond and near normal 0 2 bond would 
suggest that reactions like 

CIOO 4 X ^ CIO -f XO (6) 

should be significantly slower than similar reactions involv- 
ing H0 2 . On the other hand, the weakness of the CIO bond in 
CIOO should facilitate the extraction of the Cl atom in reac- 
tions like 

X 4 CIOO — XCl4 0 2 . (7) 

The ratio of the rates of reactions (6) and (7) when 
X = Cl as well as the comparison of the rates of the 
Cl 4* HOO and Cl 4 CIOO reactions seem to support this 
rational. Indeed, the present results suggest that CIOO 
might be best viewed as a Cl atom “chaperoned” by an 0 2 
molecule. From this point of view reaction (7) would be a 
chaperoned Cl addition rather than an abstraction reaction. 
Since the equilibrium between Cl 4 0 2 and CIOO favors the 
free Cl atom under atmospheric conditions, and since CIOO 
is expected to react like a (chaperoned) Cl atom, the CIOO 
reaction rates would have to be at least two orders of magni- 
tude faster than equivalent Cl atom reactions in order to 
compete with the latter. This might occur in cases where the 
chaperone caused a chlorine addition to be biomolecular 
rather than termolecular. 

Although there are a number of possible reactions of 
this type involving OH, H0 2 , NO, N0 2 , and N0 3 which 
might have some significance under the conditions of labora- 
tory investigations, none of them would seem to be impor- 
tant under average daylight atmospheric conditions. For ex- 
ample, the reaction between CIOO and OH yielding 
HOC1 4 0 2 would be expected to have a rate similar to that 
for the Cl 4 CIOO — ► Cl 2 4 0 2 reaction (/c — lx 10” 10 
cm 3 molecule -1 s -1 ), but the formation of HOC1 by this 
mechanism would be orders of magnitude lower than by the 
CIO 4 H0 2 reaction and the rate of removal of free Cl would 
be insignificant. The CIOO reactions with NO and N0 2 may 
compete with the third order Cl atom addition reactions, but 
in this case the products are photolyzed rapidly and are not 
thought to be important in the atmosphere. 
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Full valence and first-order Cl wave functions are invariant with respect to unitary transformations among 
the valence orbitals. We exploit this degree of freedom and show that by transforming the valence orbitals 
into a corresponding orbital basis, nonorthogonai transition moment calculations become an easily managed 
task. Sample full valence calculations on several states of 02 and OF are also presented. 


INTRODUCTION 

The computational difficulties involved in calculating 
the transition moments between Cl wave functions which 
are composed of nonorthogonai orbitals are well 
known . 1,2 These difficulties arise from the need to 
compute the cofactors of the overlap matrices associ- 
ated with most pairs of determinants comprising the 
two Cl wave functions. This problem vanishes if the 
same molecular orbitals can be employed in both Cl 
wave functions. Unfortunately, the low-lying states of 
a large number of molecules are not well described by 
the same orbital basis, e. g. , BeO and MgO . 3 In the 
case of a limited Cl expansion the nonorthogonai transi- 
tion moments can be obtained without excessive com- 
putational effort . 4,5 For larger systems, the bounds on 
a transition moment can be estimated by performing 
orthogonal calculations with both sets of orbitaLs. More 
rigorous bounds to the transition moment can be obtained 
from techniques developed by Weinhold 6 and extended by 
Langhoff and Chong . 7 

In this paper, we show that the computational difficul- 
ties encountered in nonorthogonai transition moment 
calculations involving full valence, first-order and to a 
lesser degree second-order Cl wave functions can be 
greatly reduced by exploiting the invariance of the Cl 
wave functions with respect to unitary transformations 
among the valence orbitals. In this algorithm the same 
core orbitals are employed in both states and the vaLence 
orbitals of the two states under investigation are trans- 
formed to a corresponding orbital basis . 9 As a result 
of the simple nature of the overlap matrix, the calcula- 
tion of a nonorthogonai transition moment becomes an 
easily managed task. 

The ability to compute nonorthogonai transition mo- 


a) NASA-NRC Postdoctoral Fellow 1978-1980. Present address: 
IBM Research Laboratory, 5600 Cottle Road, San Jose, Calif. 
95193. 

b) This report was prepared as a result of work performed by 
this author under NASA Contract No. NAS1-14472 and NAS1- 
15810 while in residence at ICASE, NASA Langley Research 
Center, Hampton, Va. 23665. Present address: Polyatomics 
Research Inc. 1101 San Antonio Road, Mt. View, Calif. 94043. 


ments for a complete active space (CAS) wave function 
with littLe work in excess of that needed to obtain the 
orthogonal transition moment is especially useful in 
light of the recent advances by Roos 9 and Siegbahn 10 on 
efficient methods to perform large CAS MCSCF’s. 
Moreover, recent studies of the alkaline earth oxides 
have indicated that the orbitals needed to compactly de- 
scribe the ground state and the low-lying excited states 
are significantly different . 2,3 Clearly, any method which 
can be efficiently employed to obtain the nonorthogonai 
transition moments for such systems are very useful in- 
deed. 

Corresponding orbitals have been employed extensive- 
ly in connection with projected-unrestricted-Hartree- 
Fock calculations , 11 but they have rarely been employed 
in connection with Cl calculations. The notable excep- 
tions are the work of Davidson and Martin 12 in which 
corresponding orbitals were exploited for interpretative 
purposes and the work of Bagus 13 where corresponding 
orbitals were used in the study of localized solutions in 
core ionized systems. 

King 14 has used corresponding orbitals as an alterna- 
tive to the method of Prosser and Hagstrom 15 for com- 
puting the matrix elements in nonorthogonai Cl and tran- 
sition moment calculations. In King’s method a differ- 
ent set of corresponding orbitals are determined for 
most determinant-determinant interactions. 

CALCULATION OF TRANSITION MOMENTS 

Our goal is to compute the transition moment (TM) 
between two states t/t i and tp 2 : 

TM = (il>i\x f y, or z | i/>2> . 

The most direct way to proceed is to compute the tran- 
sition density matrix 18 (TDM), and then take the trace 
of the dipole integrals and the TDM to obtain the TM. 

In practice, ip is a linear combination of CSF’s 

*/•£ ClJ.CSFf , 

j 

and each CSF is a linear combination of determinants 

CSF'=2Zc'(fc)«DET'(fc) . 
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For wave functions of this type, the contribution of each 
determinant -determinant interaction [weighted by the 
appropriate Cl/s and C^’s] is determined and added 
to the TDM. If we let m denote an occupied orbital in 
the one determinant and n denote an occupied orbital in 
the other, then each contribution to the TDM can be ex- 
pressed as follows: 

tdm( w ,k) = Eci'ci' £c'(o)C?(£)A* , 

i»j o»p 

where is the cofactor of the o,p determinant-deter- 
minant overlap matrix. 

In the most general case of nonorthogonal orbitals, 

A mn is nonzero and the contributions from each pair of 
determinants must be evaluated. The inclusion of sym- 
metry gives rise to some savings, but the amount of 
work is still very large. 

The use of orthogonal orbitals represents a large re- 
duction in work, since is either zero, one, or minus 
one, and most of the zero detection can be accomplished 
at the occupation level. Only identical occupations or 
occupations which differ by one orbital need be pro- 
cessed. Because a cofactor need not be evaluated for 
every determinant -determinant interaction, an algo- 
rithm can be designed to evaluate orthogonal transition 
moments in which very little work is performed at the 
determinantal level. 

A third possibility is to employ corresponding orbitals 
as the one-particle basis in the Cl calculations. Cor- 
responding orbitals are defined by the following overlap 
properties 8,14 : 

<a i |a j >=(6 i |6 j ) = 5 u (la) 

and 

k i b,) =<f j 6 iy . (lb) 

Given two sets of orbitals (X*) and (Y^) the correspond- 
ing orbitals {aJ associated with (X t ) are obtained by 
diagonalizing: 

G = S xy (S xy )\ (2) 

where =<X i l Yj ) and i , j are valence (active) orbitals. 
The corresponding orbitals associated with (T^) can 
be obtained by diagonalizing: 5 yr (S yx ) T . However, this 
procedure can fail to yield a diagonal overlap matrix if 
degenerated eigenvalues are present. This problem 
can be solved by diagonalizing G to obtain one set of 
corresponding orbitals and then obtain the second set 
by means of a projection operator 17 : 

£*k> I !*><** k> (3a) 

; 

= Ek>kk>=rfik> • (aw 

j 

The diagonal overlap matrix associated with the cor- 
responding orbital basis can be used to detect zero con- 
tributions to the TDM at the occupation level and also 
at the determinantal level. Moreover, the simple struc- 
ture of the overlap matrix greatly reduces the computa- 
tional effort required to evaluate the cofactors. Thus, 
in the Cl studies where corresponding orbitals can be 


employed, it is possible to perform a nonorthogonal 
TM calculation with only slightly more work than the 
effort required to obtain the comparable orthogonal 
TM. 

We should note that if one is interested in calculating 
the TM between several states, the corresponding or- 
bitals between each pair of states must be determined 
and the Cl calculations repeated in this basis. How- 
ever, we find with our Cl that the time to evaluate the 
wave functions is usually less than the time to compute 
the nonorthogonal TDM. 

FULL VALENCE Cl 

In a full valence Cl calculation, the energy is invariant 
to unitary transformation among the valence orbitals. 

One is then free to choose the orbital basis which re- 
duces the computational effort needed to obtain the tran- 
sition moment. This task can be accomplished very 
easily if the overlap matrix between the occupied or- 
bitals of the two states is diagonal. The orbitals which 
satisfy this overlap criteria are simply the correspond- 
ing orbitals. Because of the diagonal nature of the over- 
lap matrix, two occupations which differ by two or more 
orbitals do not contribute to the TDM and are not pro- 
cessed. In fact, these calculations are performed in 
the same manner as an orthogonal transition moment 
computation, the only difference being that one multi- 
plies each contribution of the TDM which would arise 
in an orthogonal calculation by a cofactor of the overlap 
matrix. These cofactors are particularly simple to ob- 
tain, as every nonzero cofactor is of the form 

CF< t) = <4a) 

= n k , (4b) 

where the product extends over all of the spin orbitals 
in one of the interacting occupations. It is convenient 
to precompute the terms in brackets in Eq. (4a) for 
every occupation appearing in one of the Cl wave func- 
tions. If one of the ^’s appearing in Eq. (4a) is zero, 
then one precomputes the product defined by Eq. (4b). 

A key is also constructed to indicate which orbital pos- 
sesses a zero overlap, as only this orbital is involved 
in the nonzero contributions to the TDM. If two or more 
of the di’s are zero, then this occupation cannot con- 
tribute to the TDM. 

Thus, in the full valence case the nonorthogonal TM 
is evaluated in virtually the same time as the orthogonal 
TM. This is very useful for CAS MCSCF wave func- 
tions. 

FIRST-ORDER Cl 

Corresponding orbitals may also be utilized to great 
advantage in the calculation of transition moments in- 
volving first-order Cl (FOCI) 10 wave functions. The 
formulas for the cofactors which occur in a FOCI cal- 
culation are naturally more complicated than those en- 
countered in a full valence calculation. However, these 
formulas represent a significant reduction of work when 
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compared to the effort required to perform the same 
calculation without employing corresponding orbitals. 

A FOCI is simply described as all single excitations 
away from a full valence Cl. The wave function is there- 
fore invariant to unitary transformations among the va- 
lence orbitals. The valence orbitals comprising the Cl 
wave function are transformed to the corresponding or- 
bital basis. The overlap between the valence orbitals 
is now diagonal, while the overlaps valence(l)-virtual(2), 
valence(2)-virtual(l), and virtual(l)-virtuai(2) are non- 
zero. The wave function is also invariant to unitary 
transformations among the virtual orbitals; however, 
a corresponding orbital transformation does not signifi- 
cantly simplify the calculation of the TDM and is there- 
fore ignored. 

From the definition of a FOCI it is clear that each 
CSF can contain at most one virtual orbital. The con- 
tribution to the TDM arising from the interaction of two 
CSF's containing only valence orbitals are determined 
by the simple procedure described above. In this dis- 
cussion, the two states under consideration do not pos- 
sess the same symmetry properties; however, the for- 
mulas necessary to calculate the TDM between two states 
of the same symmetry are a trivial extension of the six 
cases presented here. 

In the discussions which follow, “active” is used to 
describe a symmetry in which the number of occupied 
orbitals in the two CSF's differ, while symmetries 
which contain the same number of occupied orbitals are 
called “passive.” The passive symmetries contribute 
only a simple multiplicative factor to the nonzero cofac- 
tors. Also, at the determinant level the overlap matrix 


can be factored into alpha and beta components. Both 
cases are treated in the same manner and thus our ex- 
amples deal with only one set of spin orbitals. 

Our general approach is to compare each pair of oc- 
cupations in the two wave functions. If the number of 
occupied orbitals in symmetry; of occupation 1 is^N],- 
then if \N) -N ) I ^ 2 for any symmetry j, or 

£ |nj-n*I>2 - 

i 

the contribution to the TDM is zero (this is the same 
test used to detect zeros in standard nonorthogonal TM 
calculations). 

At this point the active and passive symmetries are 
identified. If a passive symmetry contains only valence 
orbitals, then, as a result of the diagonal overlap ma- 
trix, the orbitals in this symmetry must be coincident 
(a { coincident with £>*) in the two occupations or the con- 
tribution to the TDM is zero. The determinants involved 
in the two occupations are then compared. The two tests 
that were made at the occupation level to detect zeros 
are then repeated at the determinant level. For inter- 
acting determinants, the noncoincident orbitals are 
analyzed and one of the six different cases is used to 
evaluate the cofactors and obtain the contributions the 
TDM. The six cases are described in detail below. In 
order to reduce the complexity of the following formulas 
we have elminated the overall phase factors. 

Case I 

A configuration with a virtual orbital in an active 
symmetry interacts with a full valence configuration: 


Sym 1 


Sym 2 Sym 3 


Symmetry 1 


Symmetry 2 


Symmetry 3 



Symmetries 1 and 2 are active symmetries, while 
symmetry 3 is a passive symmetry. The nonzero co- 
factors of this matrix arise from the exclusion of row 
6 and one of the first four columns. These cofactors can 
be written down by inspection: 

D(zo >W i)=(ndi)(n^')(n^") > <s> 

d(zo, a [) = Cj (n^)(n<)( n *"') =d{z °> % - 


where ZO denotes the valence orbital in configuration 1 
which does not overlap with any of the orbitals in con- 
figuration 2, a { and b { denote valence orbitals, v* de- 
notes a virtual orbital in the primed symmetry, and the 
c 4 * s are overlap elements between the occupied orbitals 
of one configuration and the virtual orbitals in the other 
configuration. 

Case II 

Both the interacting configurations contain a virtual 
orbital in the same active symmetry: 


j 
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t 

a u 


Once again, the nonzero cofactors of the overlap ma- 
trix arise from the exclusion of row 6 and one of the 
first four columns. These cofactors are also very sim- 
ple to evaluate: 


D(ZO , 

K)=(n^)<(n<) - 

(7) 

D(ZO, 


(8a) 


.[/(n.-;)-E(.n<)*q(n<), 




(8b) 

D{ZO y 

a\)=( n • 

\ i*au »a‘j / \ k / 

(9) 

As \ 

ve noted above, FOCI wave functions are also 

in- 

variant to unitary transformations among the virtual or- 


bitals; however, there seems little to be gained by 
transforming the virtual orbitals to the corresponding 
orbital basis since the overlap between the virtual or- 
bitals only enters Eq. (8) in a straightforward manner. 

It is also important to note that these formulas apply to 
the case where there are one or more rows and columns 
of zero in the valence-valence portion of the overlap 
matrix. This would occur if one or more of the d { ’ s 
were zero or if there was an orbital noncoincidence in 
the valence space. 


Case III 

Two configurations interact with virtual orbitals in 
different symmetries: 



a l 



v[ 

/ 1 
a x 

/ / 
a z 

b[ 


0 

0 

c[ 

0 

0 \ 

b'z 

1 0 

di 

0 

c' z 

0 

0 \ 

b'z 

0 

0 


c's 

0 

0 

b\' 

0 

0 

0 

0 

d[' 

0 

b'z 

i 0 

0 

0 

0 

0 

d'z 

Ml' 

\ 0 

0 

0 

0 

r " 

C 1 

Cz 1 


The nonzero cofactors of the overlap matrix involve 
rows 3-6 and columns 1-4: 



(10) 


(lla) 


(lib) 

D(b\', K)=(n^) f(n<) 

(12a) 

■(WlH . 

(12b) 


(13a) 


(13b) 


Case IV 

One of the interacting configurations contains a virtual 
orbital in a passive symmetry while the other configura- 
tion contains only valence orbitals: 



t 


ZO x 

The only nonzero cofactor is 

mzo z , 2o l )=^n^)(n<)<'(n<') • 
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TABLE I. Summary* of transition moment and excitation energies for Oj. 




Nonorthogonal 

TM 

Nonorthogonal 

Orthogonal TM 

Orthogonal 

Tanaka b and 
Yoshimine 
Orthogonal TM 

( 4 2;ui 4 i;> 

R =2.2 
R= 2.4 

0. 568 
0.314 

0.2241 

0.2574 

0. 5594 
0.3109 

0. 2203 
0.2547 

0. 5536 
0. 3091 

<Xi*i 4 z;> 

R =2.4 

0. 1983 

0.0755 

0. 1920 

0.0737 

0. 1918 


*A11 quantities are in atomic units. 

^nly the transition moment results were presented in Ref. 21. This orthogonal calculation employed 4 I" orbitals. 


a\ a' z a' 3 v\ v‘ z a\* a z a’" a z " 



=S 


-zo 


The nonzero cofactors arise from excluding row ZO 
together with an orbital in symmetry 1 (i. e. , a \ , a Zi yj). 
One is then left with the problem of evaluating the cofac- 
tors arising from the active symmetry block in which the 
virtual orbitals are present. There is no advantage in 
prototyping this symmetry block and the cofactors can 
be most economically evaluated by employing the stan- 
dard biorthogonalization scheme. Once again, rapid 
zero detection is afforded in this scheme due to the sim- 
ple structure of the overlap matrix in every symmetry 
but the first. 

The computational advantages of corresponding or- 
bitals may also be utilized in an approximate second- 
order Cl scheme where one generates a set of corre- 
sponding orbitals, performs a full valence Cl calcula- 
tion to obtain a list of dominant configurations, and then 
performs a single and double excitation Cl from this 
list. This technique might be used at one point on a 
potential energy surface to compare the transition mo- 
ment obtained in a first-order calculation with the single 
and double excitation Cl. Because one is employing an 


overlap criterion to generate the list of dominant con- 
figurations, this type of Cl probably could not be used 
to generate a potential energy surface. A problem 
would be expected to arise where the two potential 
curves do not parallel one another. 

SAMPLE CALCULATIONS 

Two sets of sample calculations were run using 
MOLECULE - NONAME . 19 In addition to showing the effi- 
ciency of this method, they show that the overlap of the 
corresponding orbitals can be used as a measure of the 
importance of using nonorthogonal orbitals. 

o; 

Full valence MCSCF calculations were performed on 
the low-lying states of OJ. van Duijneveldt’s 20 lls6p 
primitive Gaussian basis was contracted (611111/3111) 
to 6s4/>. This basis was augmented with two d functions 
with exponents 2. 21 and 0. 66. The valence space con- 
sisted of the 2<j g , 2 cr u , 3cr,, 3cr u , ln u , and In, orbitals 


TABLE II. Corresponding orbital overlaps for Oj. 





2a, 


3a, 



l7r r 

< 4 z;i 4 z;> 

R = 2. 

4 

0.999 993 

0. 999 947 

0.999 583 

0. 998 691 

0.999 995 

0.999 960 

< 4 s;i 4 s-> 

R= 2. 

2 

0.999 997 

0. 999 905 

0. 999629 

0. 999 341 

0.999 991 

0.999 990 

( 4 nj 4 z;> 

R = 2. 

4 

0.999 998 

0. 999 976 

0.999 507 

0.999 071 

0.999 917 

0.999 851 


L 
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TABLE III. Full valence MCSCF and Cl results for OF. 1 TABLE IV. Corresponding orbital overlaps for OF. 


orbitals 


Ml orbitals 


z n 

Orthogonal TM 


- 174.004 317 
— 174. 240 015 
0.241098 
0. 029 800 


- 173. 953 750 
-174.257 807 
0.304 057 
0. 051 000 


Optimum =0. 25349 nonorthogonal optimum TM =0. 0334 


a All quantities are in atomic units. 


and the lo orbitals were frozen in all calculations. 

These orbitals were obtained from 6-state ( 4 TJ)SCF 
calculations. In Table I we compare our results to 
those obtained by Tanaka and Yoshimine 21 -in which they 
employed a basis of Slater orbitals. In Table II, the 
overlaps between the corresponding orbitals obtained in 
the three calculations are presented. 

The nonorthogonal and orthogonal transition moments 
are in close agreement with the results of Tanaka and 
Yoshimine. At this level of approximation, the c -state 
( 4 N~) orbitals offer a very good representation of the 
three states investigated. The very high overlaps be- 
tween the sets of corresponding orbitals show that the 
valence space is very similar for all states, so it is not 
surprising that orthogonal orbitals work well in this 
case. 

OF 

Full valence MCSCF (including the 2s orbitals) calcu- 
lations were performed on the 2 n and 2 A states of OF. 
van Duijneveldt’s llsQp Gaussian basis was used and 
contracted as before. This basis was also augmented 
by two d functions on each center with exponents 2.901, 
0.6194, 2.214, and 0. 655 on oxygen and fluorine. The 
results of the OF calculations are presented in Table 
III and the overlaps of the corresponding orbitals are 
presented in Table IV. These calculations were carried 
out in C 2v symmetry which led to Cl expansions of 42 
and 18 CSF’s for the 2 a and the 2 H states, respectively. 
All of the transition moment calculations (which were 
undertaken in the corresponding orbital basis) required 
approximately an order of magnitude less time than the 
equivalent general nonorthogonal transition moment cal- 
culation. 

In the OF calculations, the Prosser-Hagstrom non- 
orthogonal transition moment calculation required 25 
cp sec while the equivalent corresponding orbital calcu- 
lation required 2 cp sec. The orthogonal transition mo- 
ment calculation required approximately 2 cp sec as 
well. (These calculations were performed on a CDC 
Cyber 173 computer.) 

A nonorthogonal treatment of OF clearly indicated 
that the transition moment and the energy separation 
fluctuate significantly when different sets of orthogonal 
orbitals are used in the calculation. The excited state 
orbitals clearly provide a better description of both 
states than do the ground state orbitals. The overlap 


jO 4(7 jO <ia _ l - J- 

U.99999H «. 999990 0.9997*9 0.973 1^0 ('.9999*0 n. 999 T'Ct 


between the corresponding orbital reflects significant 
orbital reorganization in the excited state. Therefore, 
it is not surprising to find that the TM obtained in or- 
thogonal and nonorthogonal calculations do not coincide. 
This case is more representative than OJ where the or- 
bitals are virtually the same for both states. This ef- 
fect is expected to become even more pronounced in 
heteronuciear diatomics which possess a large ionic 
character. 2,3 

CONCLUSIONS 

Corresponding orbitals provide a very convenient basis 
in which to calculate transition moments from full va- 
lence and first-order Cl wave functions. Moreover, 
corresponding orbitals have been shown to reduce the 
work needed to compute the cofactors of the overlap ma- 
trices and to detect zero contributions to the transition 
density matrix. In the case of the first-order Cl, the 
higher the symmetry, the more efficient is the zero de- 
tection. 

This technique provides an attractive compromise be- 
tween the standard orthogonal transition moment algo- 
rithm and the nonorthogonal transition moment method 
of Prosser and Hagstrom. 15 This is especially true in 
large Cl calculations in which the states of interest are 
poorly described by the same orbital basis. 

It is also noted that the overlap between the corre- 
sponding orbitals can be used to identify the systems in 
which a nonorthogonal calculation is required. 
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